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ABSTRACT 


Results  of  three  calculations  of  the  impulsive 

interaction  of  thermal  radiation  from  high- alt i tu&e 

nuclear  explosions  with  targets  within  the  fireball 

(the  thermomechanical  effect)  are  reported.  These 

problems  constitute  a  preliminary  investigation  into 

techniques  of  calculation  and  phenomena  of  importance 

to  be  followed  later  by  a  survey  to  establish  lethal 

radius  envelopes.  The  interaction  with  200-KT  yield 

at  60,000-ft  altitude  is  found  to  be  predominantly 

the  effect  of  the  blast  wave  which  is  perhaps  just 

4 

lethal  at  a  range  of  2-5  x  10  cm.  The  interaction 

of  the  same  yield  at  9°,000  ft,  however,  does  show 

important  thermomechanical  as  well  as  blast  effects. 

The  lethal  radius  is  found  to  be  greater  than 
4 

3.0  x  10  cm.  Presentation  of  the  calculated  results 
is  supplemented  by  discussion  of  the  method  of 
calculation,  preparation  of  initial  data,  and  evalua¬ 
tion  of  assumptions  constituting  the  interaction 
model. 
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INTRODUCTION 


Preliminary  results  of  calculations  of  the  interaction  of  thermal 
radiation  from  a  nuclear  explosion  with  an  object  within  the  fireball  are 
presented.  The  effect,  called  the  thermomechanical  effect  (TME),  results 
from  the  irradiation  of  a  structure  by  intense  thermal  energy  emitted  by  the 
heated  air  and  the  subsequent  vaporization  and  blowoff  of  the  heated  vapor. 
The  reaction  of  the  blowoff  to  exert  an  impulse  on  the  structure  produces 
the  TME. 

In  this  report  three  calculations  are  described  in  some  detail.  Two  of 
these  use  the  environment  created  by  a  200-KT  yield  at  90> OOO-ft  altitude 
and  differ  in  the  range  of  the  object  (called  for  convenience  a  re-entry 
vehicle,  or  R.V.)  from  the  origin  of  the  explosion.  These  ranges  are: 

Ij,  if. 

problem  1,  2.5  x  10  cm,  and  problem  2,  3*0  x  10  cm.  The  third  problem 

4 

represents  interaction  with  an  R.V.  at  range  2.5  x  10  cm  from  a  200-KT 
explosion  at  60,000  ft. 

Material  in  the  report  is  organized  into  three  sections.  In  Section  I 
a  current  statement  of  characteristics  of  the  General  Atomic  interaction 
code,  SUPER  SHJTTER,  is  presented.  This  code  has  been  used  to  perform  the 
TME  calculations  presented  here  as  well  as  to  investigate  other  interaction 
phenomena.  Section  II  contains  an  analysis  of  the  explosion  environment  to 
obtain  boundary  conditions  for  the  TME  calculations.  A  critique  is  also 
presented  of  some  of  the  assumptions  made  about  the  environment.  Finally, 
in  Section  III  the  calculations  are  analyzed  in  some  detail  to  present  the 
characteristics  of  the  TME  interaction  process. 

At  60,000-ft  altitude  it  is  found  that  the  effect  of  ambient  fireball 
pressure  is  much  more  important  than  ablation  pressure  in  determining  the 
impulsive  loading.  Consequently,  this  altitude  is  below  the  altitude  where 
the  TME  plays  a  role  and  is  properly  considered  in  the  blast  regime .  At 
90,000  ft  both  ambient  and  ablation  pressures  are  important.  This  is  a 
transition  region  between  effects.  The  lethal  range  appears  to  be  greater 
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than  3  x  10  cm.  An  interesting  and  important  effect  of  abrupt  heating  and 
ionization  of  the  R.V.  vapor  is  identified  in  the  problems  at  90,000  ft. 

Its  quantitative  description  may  depend  sensitively  on  low-temperature 
radiative  properties. 

While  calculations  presented  here  by  no  means  constitute  an  exhaustive 
survey  of  the  TME,  it  is  thought  that  presentation  of  early  results  will 
accomplish  the  following  results: 

1.  Acquaint  a  wider  audience  with  techniques  available  for  application 
to  the  TME  and  similar  problems. 

2.  Document  more  currently  and  in  greater  detail  the  calculational 
methods  and  material  properties  in  use. 

3.  Prompt  criticism  of  the  approach  to  the  calculation  of  the  TME. 

Additional  calculations  are  now  planned  to  detail  the  lethal  radius  as  a 
function  of  altitude.  Further  improvements  in  code  and  environment  specifi¬ 
cation  are  also  being  incorporated. 
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I.  THE  SUPER  SPUTTER  CODE 

Calculations  of  the  TME  have  used  SUPER  SPUTTER,  a  general  purpose,  one¬ 
dimensional  radiation  transport  hydrodynamics  code  baaed  on  the  Lagrangian 
description  of  the  fluid.  Several  modifications  of  this  code  have  been  made 
to  carry  out  ‘lME  calculations.  In  the  following  sections  the  most  important 
of  these  code  changes  are  described  in  order  to  illuminate  the  nature  of  the 
physical  model  and  some  of  the  mathematical  approximations  being  used. 

In  section  A  the  photon  source  treatment  is  described  in  complete  detail; 
section  B  indicates  how  the  ambient  pressure  inside  the  fireball  is  applied 
as  a  boundary  condition  to  the  vapor;  section  C  describes  the  treatment  of 
transition  between  vapor  and  solid  at  the  "boiling"  interface  and  the 
divergence  of  expanding  vapor;  section  D  indicates  how  deviations  from  local 
t hemodynamic  equilibrium  which  may  occur  in  optically  thin  and  low-density 
regions  are  taken  into  account;  section  E  describes  the  radiation  transport 
routine  used  in  the  code;  section  F  gives  a  brief  description  of  newly 
incorporated  features  of  the  problem  generator  code  which,  in  addition  to 
preparing  initial  data,  allows  redistribution  of  zones  during  the  course  of 
the  problem.  Finally,  section  G  describes  the  properties  assumed  for  the 
target  material  and  presents  previously  unpublished  calculations  for  some 
of  the  temperature -density  regime. 

The  successful  completion  of  a  number  of  problems  described  below  has 

both  tested  the  modifications  mentioned  herein  as  well  as  shown  up  some 

inconsistency  in  the  angular  distribution  used  in  an  earlier  calculation •  ^ 

The  new  routines  have  been  thoroughly  checked  and  are  producing  consistent 

results  with  obvious  improvements.  In  one  case  a  comparison  of  two  runs 

differing  in  the  angular  distribution  of  the  source  which  resulted  from  an 

earlier  incorrect  source  formulation  is  presented.  The  error,  contained  in 
(?) 

our  first  report v  '  on  the  thermomechanical  source,  has  been  corrected  in 
the  description  of  the  first  section  and  does  not  occur  in  the  remainder  of 
problem  results  presented. 

A  more  complete  report  on  the  SUPER  SPUTTER  code  is  now  in  preparation, 
with  publication  expected  shortly. 
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A*  THERMAL  SOURCE 

The  SUPER  SPUTTER  code  will  he  used  to  calculate  the  interaction  of  a 
high-temperature  radiative  source  on  an  object  placed  within  the  thermal 
environment  o  The  description  of  the  problem  is  most  closely  approximated  by 
that  for  the  single -frequency  source  used  in  the  investigation  of  optical 
interactions « ^ 

A  generalization  of  the  optical  interaction  source  has  been  made  to 
describe  the  thermal  interaction.  In  it  the  incident  radiation  is  assumed 
to  fall  on  the  locally  plane  surface  of  the  target  or  target  vapor.  The 
source  is  characterized  by  a  total  intensity  which  may  change  with  time,  a 
frequency  distribution  approximated  by  representative  groups  of  photons,  and 
one  of  several  different  choices  of  angular  distribution.  The  photon  groups 
are  absorbed  in  the  target  vapor  according  to  an  approximation  to  the  ionic 
absorption  coefficient  which  takes  into  account  free-free  and  bound-free 
transitions  of  the  ionized  vapor.  This  coefficient  depends  on  the  photon 
frequency  u,  the  material  temperature  9,  and  density  p,  as  well  as  on  the 
target  material  and  degree  of  ionization.  The  absorption  of  each  frequency 
group  in  the  vapor  is  determined  by  numerical  integration  to  find  the  optical 
depth  from  the  surface  corrected  for  induced  emission  at  the  frequency  in 
question.  The  flux  is  then  obtained  as  a  function  of  the  optical  depth  and 
parameters  specifying  the  nature  of  the  incident  angular  distribution . 

The  Source 

The  thermal  source  is  specified  by  the  angular  and  frequency  distributions 
of  the  photons  incident  on  the  target.  While  in  the  most  general  problem 
these  distributions  cannot  be  factored  and  depend  on  the  time,  in  this  first 
approach  several  simplifying  assumptions  have  been  made.  Perhaps  most 
fundamental  of  these  is  that  the  frequency  and  angular  distributions  may  be 
factored  into  a  product  of  two  individual  distributions,  one  of  which  depends 
only  on  frequency  and  the  other  only  cn  angle. 

Angular  Distribution 

The  source  routine  now  offers  three  options  for  the  photon  angular 
distribution,  each  of  which  is  a  simple,  time -independent  function  for  which 
angular  groups  are  unnecessary.  The  choices  are i 
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!•  Isotropic  distribution,  l(Q)  =  I  . 

2.  Cos  distribution,  l(0)  =  I 

3*  Ray  or  pencil  distribution,  l(Q)  =  Io8(Q  •  tJ0) 

In  the  second  distribution  }i  is  the  cosine  of  the  angle  from  the  ray 
direction  to  the  normal  to  the  assumed  plane  target.  In  the  last  distribution 
a  pencil  of  photons  falls  on  the  target  in  the  Dq  direction  for  which  m-  = 

All  of  these  distributions  allow  the  planarity  of  the  interaction  to  be 
preserved  and  lead  to  the  following  simple  functions  f  (x)  describing  the  net 
flux,  0  =  f(ft)  cos  9  dO,  of  the  source  photons  at  optical  depth  x  normal  to 
surface  within  the  target: 

1.  f  (x)  =  2E3(x) 

2.  f(x)  =  3E^(x) 

3.  f(x)  =  exp  (-x/u0) 

Each  of  the  distributions  is  normalized  to  1  at  x  =  0  if  the  distribution  is 
symmetrical  about  the  normal.  For  case  3  and  oblique  incidence  (hq  ^  1), 
the  coefficient  takes  into  account  the  increased  target  area  intercepted 
by  a  beam  having  unit  area  normal  to  the  ray. 

The  first  option  is  relevant  to  the  thermal  source  inside  an  optically 
thick  fireball  where  isotropy  of  the  radiation  obtains .  Well  outside  the 
fireball  radiation  comes  from  a  rather  localized  direction  represented  by 
the  pencil  of  option  3«  Option  2  is  an  intermediate  case  perhaps  appropriate 
near  the  edge  of  the  fireball. 

Addition  of  further  simple  options  is  easily  made.  Angular  grouping  for 
a  complex  distribution,  however,  would  present  a  harder  but  still  tractable 
job. 


Frequency  Distribution 

Photons  in  the  Incident  beam  are  described  by  a  distribution  function 
f(u)  which  is  assumed  to  depend  on  time  but  not  on  angle.  The  source  calcula¬ 
tion  approximates  this  distribution  in  the  frequency  group  approximation 
wherein  the  photons  falling  within  the  interval  (u  ,  un+^)  according  to  the 
absolution  coefficient  of  a  representative  frequency.  In  the  general  case 
the  intensities  An  within  these  groups  would  be  permitted  to  depend  on  time 
in  a  quite  arbitrary  way  as,  for  example,  when  specified  by  piece -wise  linear 
expressions . 
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In  the  special  case  of  the  black  body  function,  f(u)  «  \z?‘  ■ —  ,  the 

ioo  ^ 

frequencies  may  be  chosen  in  such  a  way  that  considerable  JQ  v 
simplification  results  through  maintaining  a  constant  time- independent  ratio 
between  groups  even  though  the  black  body  source  temperature  9q  may  change 
with  time.  This  Is  accomplished  by  the  assignment  un  «  unGo^  in  ^e 

ur  are  constants  chosen  to  give  the  desired  accuracy  in  the  source  description. 
The  frequency  group  boundaries,  consequently,  depend  on  time  through  time- 
dependence  of  0^.  A  representative  choice  of  frequency  groupB  and  weights 
for  a  black  body  is  given  in  Table  1.  The  are  normalised  to  give  unit 
sum. 


Table  1.  Nine-group  Constants 

Group  Boundary  Average  Energy 

n  hu/9  An  E  =*  hu/9 

o  n  o 


0 


1 

1.0 

0.05 

0.5 

2 

2.0 

0.14 

1-5 

3 

3.0 

0.21 

2.5 

k 

4.0 

0.20 

3-5 

5 

5.0 

0.16 

4.5 

6 

6.0 

0.10 

5-5 

7 

8.0 

0.10 

7.0 

8 

10.0 

0.03 

9.0 

9 

12.0 

0.01 

11.0 

In  the  currently  available  version  of  the  source  the  group  constants 

are  loaded  into  the  code  from  Table  1  for  application  to  the  black  body 

4 

source.  The  black  body  intensity,  ,  is  allowed  to  change  with  time 
through  specification  of  constants  in  a  piece -wise  linear  variation  of  9^ 
with  time.  Generalization  of  the  routine  to  calculate  non-black  body  cases 
is  made  trivially  if  the  intensity  Is  independent  of  time  and  requires 
additional  coding  to  provide  the  description  of  the  time  dependence  of  the 
intensities  for  the  more  general  problem. 
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Source  Energy  Deposition 

The  amount  of  energy  q^  [ergs/(cm2sec)  ]  deposited  in  each  of  the  zones 
of  the  vapor  is  hounded  by  Lagrangian  coordinates  x^  and  x^+^.  The  source 
is  obtained  by  taking  the  difference  between  the  net  fluxes  incident  on 
and  emergent  from  the  zone  in  question. 

4i  ^i+1 


The  flux  is  formed  from  the  finite  difference  approximation  to  the  following 
expression  for  the  total  net  flux  integrated  over  frequency  and  angles: 


0(x) 


where  E  is  the  integrated  intensity  falling  on  the  target  and  x^  is  the 

optical  depth  at  frequency  u  measured  normal  to  the  target  surface 
x 

Xu  ~  Jo  P^^x^dx1 .  In  the  finite  difference  approximation  for  the  frequency 

groups  and  zones,  this  becomes  0.  =  E  E  A  fix  ),  where  x  is  the  optical 

i  n  n  J  n  n 

depth  for  the  nth  group  and  f  can  take  the  form  of  the  three  options  given 

J 

above . 

For  the  black  body  source 

x. 

i 

01  =  2Efg  AnE3[J  Kn(x)pdx 

o 

where  the  optical  depth  is  approximated 
x.  i-1 

JoKn(x)odx  =  XpjKC'vVV  '  (XJ+1‘XJ)> 

0  j  =  l 

•  4  /  2 

E  =  a ,  and  f  is  a  constant  to  allow  for  1/r  reduction  in  intensity. 


Flow  Diagram 

The  specification  given  above  has  been  incorporated  into  the  "QUE" 
subroutine.  This  routine,  one  among  several  options  which  may  be  incorporated 
into  SHJTTER,  is  selected  for-  thermal  mechanical  effect  calculations  while 
other  subroutines  remain  largely  in  standard  form.  The  logical  flow  of  the 
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code  subroutine  is  given  in  Appendix  A  of  this  report,  together  with  the 
FORTRAN  instructions  corresponding  to  it.  Both  flow  diagram  and  code  apply 
to  the  black  body  source.  Addition  of  a  specification  of  the  way  in  which 
the  frequency  distribution  changes  with  time  will  be  required  for  the  more 
general  frequency-dependent  problem. 


Vapor  Absorption  Coefficient 

The  Raizer  approximation  to  ionic  absorption  coefficient v  *  is  used 
to  obtain  the  dependence  on  frequency,  material  temperature,  density  and 
species.  The  prescription  including  induced  emission  is 


Kff(eU-l) 

Kff(l-e"u)eUK 


if  u  <  u„ 


if  u  i  uT, 


where 


K  =  8.62  x  1010  -§  - I 

ff  A2(hv)V 


^  =  min  (uL,  inJT) 

,  ,,  1J.6  Z? 

U1  ”  K  ’  shell  energy  =  - g - 

f  =  electron  degeneracy  parameter  = 


10“2A93/2 


Other  quantities  appearing  are 

A  =  atomic  number  of  target  material 

o 

p  =  density  (gm/crrr) 

9  =  temperature  (ev) 
hu  =  phot on  energy  (ev) 

Z  =  mean  number  of  free  electrons  per  target  atom  (obtained  in  the 
equation  of  state  subroutine) 
u  =  hv/9 


zs  - 


if  Z  < 


Z  if  Z  £  1 


Mean  square  number  of  free  electrons 
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Future  Applications 

Further  generalization  and  improvement  of  the  interaction  code  would 
indicate  modification  in  at  least  two  respects: 

1.  A  more  general  angular  distribution  would  require  analysis  to 
determine  how  the  angular  integrals  may  be  performed.  It  is  likely  that  a 
simple  distribution  could  be  handled  analytically  while  a  more  complicated 
one  would  require  angular  grouping. 

2.  Optical  properties  in  the  ionic  region  exist  for  some  target 
materials  in  greater  detail  than  the  Raizer  approximation.  For  these  (e<>g., 
CH^,  air.  Be,  Si,  Al)  the  deposition  can  be  more  accurately  performed  in 
principle.  In  practice,  however, '  this  treatment  contains  many  of  the  limita¬ 
tions  of  the  many -frequency  radiation  transport,  primarily  in  lengthy,  slow 
calculation.  It  is  probable  that  improved  material  properties  can  be  used  in 
the  source  without  excessive  lengthening  of  the  code*  The  low  temperature 
material  properties  (solid  and  molecular  regimes)  are  known  very  poorly, 
constituting  a  limitation  for  applications  having  small  incident  radiant  flux. 
Additional  work  is  needed  for  this  use  as  well  as  for  the  optical  interaction 
application  at  low  incident  power. 

B.  PRESSURE  BOUNDARY  CONDITION 

Air  inside  the  fireball  heated  to  its  ambient  temperature  exerts  a 
pressure  which,  increasingly  at  lower  altitudes,  influences  the  hydrodynamics 
of  the  vapor  blowoff.  When  the  air  pressure  i~  comparable  with  or  exceeds 
the  ablation  pressure,  me,  the  blowoff  is  inhibited,  the  density  remains 
high,  and  vapor  opacity  is  affected.  This  pressure  also  increases  the 
impulse  experienced  by  the  target.  If  the  velocities  of  the  air  and  vapor 
remain  small  compared  with  thermal  velocities,  the  effect  of  the  air  is  to 
apply  a  static  pressure  to  the  air-vapor  interface,  constituting  a  boundary 
condition.  Further  discussion  of  this  point  is  found  in  section  II. D.  At 
a  particular  radial  distance  from  the  explosion  point  this  pressure,  given 
by  the  ambient  fireball  pressure,  changes  with  time  as  shown  typically  in 
Fig.  II. 2. 

Modification  of  the  SUPER  SPUTTER  code  has  been  made  in  the  HYDRO 
subroutine  to  accommodate  the  externally  applied  boundary  pressure.  The 
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"boundary  condition  at  the  vapor-air  interface  is  given  in  the  form  of  a 
piece -wise  linear  fit  in  time  to  the  ambient  fireball  temperature*  Flow 
diagram  and  code  are  presented  in  Appendix  B. 


C.  BOILING  ANT)  VAPOR  DIVERGENCE 

In  order  to  describe  the  transition  from  solid  to  vapor,  called  "boilingi* 

it  is  assumed  that  a  "boiling"  temperature,  9^,  and  a  constant  specific  heat, 

C„,  can  he  defined  for  the  solid*  In  turn  a  specific  energy  for  "boiling”, 

H_  =  C,J0L*  is  then  known.*  When  an  element  of  mass  is  heated  to  9,  ,  it  is 
C  V  b  d 

assumed  that  an  additional  amount  of  specific  energy  called  the  vaporiza¬ 
tion  energy,  is  required  before  the  vaporization  will  occur  and  the  temperature 
can  change.  For  E  >  +  H^.,  the  material  is  described  in  terms  of  vapor 

equations  of  state. 

As  vapor  is  evolved  at  the  solid-vapor  interface,  it  expands,  dropping 
in  density.  If  It  continued  to  expand  in  only  one  direction,  a  steady  state 
could,  not  be  achieved  in  the  vapor.  In  fact  the  vapor  will  tend  to  expand 
parallel  to  the  target  surface  also.  In  order  to  take  this  effect  into 
account,  a  mass  depletion  scheme  has  been  added  to  the  code. 

The  radial  divergence  can  be  approximately  taken  into  account  in  a 
Lagrangia n  plane  geometry  calculation  if  the  mass  in  a  given  zone  is  reduced 
as  the  zone  expands.  A  parameter  r,  the  radius  of  the  target  or  target  area 
in  the  case  of  a  laboratory  experiment,  is  introduced  into  the  calculation. 

The  mass  of  a  given  zone  on  each  time  cycle  is  then  said  to  he 


or 


in 


where 


jft 


—  (  n 

r  +  c(zx 


On 

) 


and  z  is  the  average  position,  measured  perpendicularly  to  the  target 
X 

surface  at  time  zeros  of  the  zone  i  at  time  n.  This  treatment  assumes  that: 
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(l)  the  vapor  diverges  parallel  to  the  target  surface  at  the  same  rate  at 
which  it  expands  perpendicularly  to  the  surface;  (2)  the  temperature  and 
density  of  the  material  are  functions  of  z  and  t  only;  and  (3)  vapor  which 
flows  out  of  a  tube  perpendicular  to  the  target  surface  and  of  radius  r  can 
he  removed  from  the  calculation  and  ignored.  The  work  done  and  the  total 
energy  removed  from  the  calculation  as  this  material  is  removed  frcm  the 
zone  is  taken  into  account. 

D.  DEVIATIONS  FROM  LOCAL  THERMODYNAMIC  EQUILIBRIUM 

In  the  very  low  dens  it  y  leading  edge  of  the  blownoff  vapor,  the  material 
becomes  optically  thin.  When  this  occurs  the  local  equilibrium  between  the 
material  and  the  radiation  field  is  maintained  only  if  the  rate  of  increase 
of  the  internal  energy  by  collisional  ionization  (Ecq^)>  exceeds  the  rate  of 
decrease  by  radiation  recombination  (E  ) .  This  breakdown  of  the  local 
equilibrium  between  the  radiation  field  and  the  material  will,  in  general, 
occur  before  the  random  kinetic  energy  of  the  ions  and  electrons  fall  out 
of  local  equilibrium  with  one  another*  If  it  is  assumed  that  the  ions  and 
electrons  are  in  equilibrium,  then  separate  rate  equations  for  the  internal 
energy  due  to  random  kinetic  energy  of  the  ions  and  electrons  (R)  and  the 
internal  energy  due  to  ionization  and  excitation  (i)  can  be  integrated  if 
the  material  is  described  in  terms  of  three  independent  variables :  tempera¬ 
ture,  density,  and  the  mean  number  of  free  electrons  per  average  ion  (Z) . 

The  correct  radiation  rate  can  be  expressed  in  terms  of  9,  p,  and  Z  also* 

In  order  to  do  this,  I  as  a  function  of  Z  for  the  material  in  question  must 
be  known.  In  the  numerical  calculation,  DTE  is  assumed  in  a  mass  element 
•until  the  material  is  optically  thin  and  E  c  ^  E  From  that  point  on, 

E  ^  is  ignored  and  separate  equations  for  R  and  i  are  integrated  in  terns 
of  p,  9,  and  Z.  The  radiation  rate  due  to  recombination  and  bremsstrahlung 
is  calculated,  and  the  absorption  of  energy  due  to  self -radiation  is  ignored. 
This  correction  in  general  reduces  the  amount  of  energy  radiated  out  of  the 
vapor. 
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E.  RADIATION  TRANSPORT 

(6 ) 

The  radiation  transport  subroutine  FLASH;  used  in  plane  geometry 
interaction  calculations,  has  been  recoded  to  make  more  accurate  and  faster 
use  of  the  exponential  integrals.  This  routine  has  also  been  coupled  with 
a  diffusion  calculation  for  optically  thick  zones  in  order  to  speed  up  the 
calculation,  in  particular,  for  the  multifrequency  time -dependent  code 
expected  to  be  operational  soon.  The  present  multifrequency  code  does  a 
"snapshot”  calculation  using  group  absorption  coefficients  stored  on 
magnetic  tape. 

The  multifrequency  calculation  brings  into  core  storage  a  single 
frequency  at  a  time.  This  calculation  is  repeated  until  all  groups  have 
been  processed.  The  calculation  will  be  performed  as  a  part  of  the  regular 
SUPER  SRJTTER  cycle  to  perform  a  truly  frequency-dependent  calculation; 
however,  in  practice  the  calculation  is  at  present  too  slow  to  do  more  than  a 
few  successive  cycles  for  a  test.  These  snapshot  calculations  are  used 
primarily  to  evaluate  radiative  properties  of  the  system  and  to  compare  with 
(and  hopefully  correct)  the  mean  absorption  coefficient  results  used  in  the 
time -dependent  calculation. 


F.  PROBLEM  GENERATOR  CODE  AND  REZONE 

1.  SUPER  SHJTTER  Variables  for  Generating  and  Rezone 


Ri 

DEHTAR. 

l 

G. 

l 

THETA. 

l 

RD. 

l 


A. 

l 


asq. 


M\ 


p,  = 


E. 

l 


El. 

l 


distance  in  centimeters  from  the  origin 

R.  , -  R. 
l+l  i 

mass  in  grains  at  i  +  J  (for  sphere  is  =  (zone  mass)/(^— ) 
temperature  in  electron  volts  at  i  +  J 
velocity  in  cm/sec  of  R 

o 

specific  volume  at  i  +  |  in  cm/gm  (plane)  or  cnr/gm  (sphere) 

2 

cross  sectional  area  at  i  in  cm 

A.2 

1 

(material  +  radiation  +  viscous)  total  pressure  at  i  +  in 

o 

ergs/cm  or  ergs/cnr 

specific  internal  energy  from  equation  of  state  at  i  +  ^ 
in  ergs/gm 

specific  internal  energy  from  integration  at  i  +  in  ergs 
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SMLR^  »  radius  of  depleted  mass 
IMDA  =  zone  of  constant  material  properties. 

2 .  Generating 

The  above  variables  are  generated  for  slab  or  spherical  geometry. 
Options  of  zoning  by  mass  or  by  radius  are  available.  In  each  case  a 
constant,  a,  is  specified  giving  the  ratio  of  either  mass  or  radius  increment 
to  the  corresponding  quantity  for  neighboring  zones.  This  quantity  is 
constant  within  a  region  but  has  possibly  different  values  in  different 
regions.  Combinations  of  the  quantities  of  a  number  of  zones,  interface  mass 
or  DELTAR,  may  be  specified  to  complete  the  problem  setup.  Among  these  the 
possibility  of  matching  G  or  DEDTAR  across  IMDA  interfaces  exists. 

The  flow  of  the  solution  is: 

(1)  Input  cards  are  read  which  describe  the  density,  temperature 
and  velocity  at  several  discrete  positions  in  the  system.  Then  the  coeffi¬ 
cients  of  a  segmented  function  fit  (constant  through  cubic)  are  calculated 
for  the  density  profile,  THETA  profile,  and  RD  profile. 

(2)  An  input  IMDA  zone  card  is  read  which  establishes  for  the  IMDA 
zone:  (a)  the  left  and  right  limits  (cm),  (b)  a  description  of  how  mass 
points  are  to  be  generated,  and  (c)  material  constants.  Then  3,  4  and  5 
below  are  executed,  finally  followed  by  editing  and  initiation  of  calculation. 

(3)  Whether  generating  or  rezoning,  given  the  left  and  right 
boundaries  (cm)  there  are  six  quantities  available  to  describe  how  the 
generated  mass  points  shall  be  related  to  each  other.  Any  two  .of  the  follow¬ 
ing  (a) -(d)  constitute  a  complete  prescription,  and  are  a  means  of  describing 
what  the  ratio  of  G±+1/0±  shall  be.  (a)  The  mass  of  the  leftmost  mass  point, 
(al)  or  the  length  of  the  leftmost  mass  point ;  (b)  the  mass  of  the  rightmost 
mass  point,  (bl)  or  the  length  of  the  rightmost  mass  point,  (c)  the  constant 
ratio  of  G  ^/g^,  (d)  the  number  of  mass  points  to  generate.  To  equate  (a) 
or  (al)  with  the  unknown  previous  last  mass  point  gives  only  one  of  the 
specifications  (b),  (bl),  (c),  or  (d) .  Special:  if  the  prescription  if 
(al),  (bl),  or  just  (bl),  the  ratio  of  the  mass  points  is  then  by  DEDTAR 
rather  than  th6ir  G.  G^,  DELTAR^,  and  R^  are  calculated. 

(4)  TBETA^  are  interpolated  from  the  THETA  profile. 
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(5)  RD^  are  interpolated  from  the  RD  profile. 

(6)  SV±,  Ai,  ASQ^  are  calculated. 

(7)  E.. ,  EI^,  are  calculated. 

(8)  SMLR^  (for  mass  depletion)  are  transferred  from  an  input 
storage  array. 

3*  Re zoning 

The  rezoning  process  first  establishes  a  step  function  density 
profile  for  the  problem  to  be  rezoned*  A  linear  segmented  fit  is  calculated 
for  the  RD  profile,  and  another  for  the  SMLR  profile  if  there  is  mass 
depletion. 

Rezone  input  cards  are  read  and  facilities  are  considered  for:  deleting 
or  adding  mass;  continuing  or  rezoning  old  mass  points;  prescription  for  how 
to  re  zone  old  mass  points  or  generate  new  mass  points;  whether  to  delete, 
add,  alter,  or  preserve  IMDA  interfaces;  whether  to  add  a  pusher.  (3)  above 
is  utilized  for  rezoning  and  adding  mass  points  and  RD  and  34LR  are  inter¬ 
polated  from  their  profiles.  Unchanged  mass  points  retain  their  original 
values  of  the  variables. 

When  rezone  input  is  exhausted,  (6)  is  employed.  Rezoned  mass  points 
have  a  calculation  of  THETA^  (conserving  old  internal  energy).  New  mass 
points  get  their  THETA^  fran  rezone  input.  E^,  EI^,  are  calculated  for 
new  and  re zoned  mass  points. 

G .  MATERIAL  PROPERTIES 

Calculations  of  the  TME  presented  herein  required  thermodynamic  and 
optical  properties  for  the  target  material  properties  for  the  target  material. 
The  properties  for  polyethylene  with  composition  (CHg)n  as  representative  of 
an  ablative  material  having  low  thermal  conductivity  and  large  specific  heat 
have  been  used.  Two  regimes  are  to  be  distinguished:  that  in  which  the 
temperature  is  low  and  ionization  is  small,  and  the  ionic  region  at  higher 
temperature . 

In  the  low  temperature  region  a  simplified  molecular  equation  of  state 
is  used.  Below  the  vaporization  temperature  of  0.125  ev  an  average  specific 
heat  of  2.4  x  1011  ergs/(gm  ev)  applies.  Thermal  conductivity  also  is 
effective  only  in  the  solid  where  the  conductivity  is  taken  to  be 
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10  ^  cal/(cm  sec  °K) .  The  heat  of  vaporization  is  the  small  value 
1o26  x  1011  ergs/gm,  sufficient  to  depolymerize  the  plastic  into  CgH^.  With 
subsequent  heating  dissociation  into  atomic  C  and  H  occurs  in  an  abrupt  but 
continuous  transition  whose  transition  temperature  depends  weakly  on  the 
density,  ^with  a  value  of  ~  0.4  ev  for  the  problems  of  interest .  A  disso¬ 
ciation  energy  of  8„5  x  lO1^"  ergs/gm  is  usedo  Optical  properties  in  the 
region  below  1  ev  are  based  on  extrapolation  from  calculations  relevant  to 
air,  not  taking  into  account  the  molecules  and  molecular  ions  actually 
present  in  CH^  at  temperatures  below  1  ev.  Correspondingly  schematic 
representations  of  the  Rosseland  and  Planck  mean  opacities  are  usedo  Further 
work  is  needed  to  remedy  this  deficiency  for  quantitative  calculations  on  CHg. 
It  should  be  emphasized  that  the  particular  material  composition  is  much 
more  critical  in  this  low  temperature  region  than  in  the  ionic  regime  at 
higher  temperature. 

In  the  ionic  region  the  thermodynamic  and  optical  properties  have  been 

(7) 

calculated  by  the  ionic  model.  Results  of  unpublished  calculations  for 

(Q) 

CHg  using  the  DIAPHANOUS  codev  ;  have  been  incorporated  into  SPUTTER* 

Taking  into  account  the  line  contribution  to  the  optical  absorption  coeffi¬ 
cient,  the  results  represent  an  accurate  estimate  of  optical  properties  of 
CH^  as  contained  in  the  Rosseland  optical  mean  free  path  used  in  the  diffusion 
regimes  and  the  spectral  absorption  coefficient  used  in  the  multifrequency 
"snapshot"  calculations  described  above*  While  the  absorption  coefficients 
have  not  yet  been  published,  the  thermodynamic  properties  and  the  mean 
optical  properties  as  calculated  in  DIAPHANOUS  are  presented  in  the  table  of 
Appendix  C. 
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II.  EXPLOSION  ENVIRONMENT  OF  THE  TARGET 

In  this  section  we  present  the  environmental  conditions  assumed  for  the 

evaluation  of  the  TME.  Calculations  of  the  effects  of  explosions  of  the 

XW50Y2^10)  at  the  altitudes  of  60,000  ft  and  90,000  ft  provide  histories 

of  envir omental  quantities  at  all  radii  of  interest.  Using  these  data, 

three  calculations  have  been  performed  for  the  following  altitudes  of 

explosion  h,and  radii  or  target  from  explosion  center,  r. 

4 

Problem  1.  h  -  90,000  ft,  r  =  2  =  5  x  10  cm. 

4 

Problem  2.  h  =  90,000  ft,  r  =  3*0  x  10  cm. 

Problem  3-  h  =  60,000  ft,  r  =  2«5  x  10^  cm. 

We  have  assumed  throughout  the  histories  of  the  interaction  that  the 
relevant  environmental  quantities  are  the  fireball  local  static  pressure  and 
isotropic  black  body  thermal  flux.  We  present  these  quantities  in  Section 
IIA  for  the  three  cases.  In  Sections  B  and  C  we  examine,  first,,  the  validity 
of  the  black  body  assumption  and,  second,  the  arguments  for  neglecting  effects 
associated  with  motion  of  air  or  target. 

A.  /mBIENT  TEN1PEKATURE  AND  PRESSURE 

In  Fig.  II-l  we  show  temperature  vs  time  and  in  Fig.  II-2  pressure  vs 
time  at  the  fixed  radial  distances  corresponding  to  the  three  cases  described 
above.  These  quantities  are  obtained  as  the  result  of  saue  interpretation 
of  quantities  interpolated  in  position  and  plotted  in  time  from  the  cited 
explosion  calculations. 

The  structure  exhibited  by  the  temperature  curves  is  interpreted  as 
representing  two  phenomena  in  the  explosion.  At  the  distances  chosen,  a 
sharply  defined  shuck  has  formed  in  the  air  and  represents  the  first  distur¬ 
bance  reaching  the  point  in  question.  Associated  with  it  is  a  rise  in 
temperature  corresponding  to  the  material  state  resulting  from  the  strong 
shock  having  a  particular  shock  speed.  While  a  slight  temperature  rise  may 
occur  in  the  material  behind  the  shock,  the  rise  is  slight  compared  with  the 
rather  abrupt  rise  associated  with  the  transit  of  the  radiative  front.  This 
disturbance  moves  forward  primarily  by  convection  of  the  shocked  air  and 
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secondarily  by  diffusion  of  the  front  into  additional  cold  material.  The 
second  rise  leads  rather  quickly  to  the  isothermal  interior  of  the  fireball 
which  is  seen  to  be  considerably  hotter  than  the  shocked  air.  Since  the 
radiative  flux  increases  as  the  fourth  power  of  temperature ,  it  is  seen  that 
the  second  rise  leads  to  fluxes  almost  two  orders  larger  than  in  the  shock. 

Modifications  of  the  explosion  calculation  data  were  found  necessary  to 
detail  the  duration  of  the  shock  plateau.  Continuous  mathematical'  rendering 
of  th^  shock  transition  by  the  artificial  viscosity  leads,  with  limited 
numbers  of  zones,  to  washing  out  of  the  shock  front.  In  each  case,  the 
shock  front  was  identified  and  replaced  in  the  curves  by  the  more  correct 
discontinuous  profile.  Duration  of  shock  temperature  was  calculated  by 
measuring  plateau  width  and  material  velocity,  a  more  accurate  way  than 
available  with  limited  printouts  in  time.  The  very  small  preheating  of  air 
ahead  of  the  shock  by  the  emitted  thermal  flux  was  neglected,  as  was  the 
flux  itself  before  shock  arrival. 

B.  RADIATIVE  ENVIRQMENT 

The  calculations^' 10 ^  of  fireball  early  history  used  to  obtain  black 
body  temperature  also  permit  more  detailed  analysis  to  evaluate  the  validity 
of  the  black  body  assumption.  At  selected  intervals  during  the  explosion 
calculation  frequency-dependent  analysis  calculations  have  been  made  deter¬ 
mining  the  radiative  flux  and  energy  density  throughout  the  fireball  interior. 
Comparisons  between  the  frequency-dependent  radiation  energy  density  and  the 
corresponding  black  body  quantity  indicate  the  degree  to  which  the  black  body 
distribution  is  justified.  In  Fig.  III-3  these  quantities  are  compared  at 
several  fireball  radii  8.3  x  10~^  sec  after  explosion  at  h  =  90,000  ft. 

While  a  complete  analysis  would  call  for  the  comparisons  at  successive  times 
at  the  same  position,  qualitatively  the  same  information  is  achieved  by 
choosing  positions  (a)  outside  the  shock  front  given  in  zone  53  corresponding 
to  time  just  before  the  shock  wave  arrives,  (b)  inside  the  shock  plateau 
given  in  zone  46  corresponding  to  the  early  shock  heated  history,  (c)  in  the 
transition  region  in  the  radiative  front  given  in  zone  39  corresponding  to 
the  temperature  rise  from  convection  of  the  radiative  front,  and  (d)  in  the 
isothermal  sphere  plateau  given  by  zone  33  corresponding  to  the  post  tempera¬ 
ture  maximum  period.  One  may  conclude  from  Fig.  I II -3  that  in  general  the 
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■black  body  assumption  is  well  justified.  In  the  plateau  regions,  indeed, 
it  is  highly  accurate,  wfrUe  in  the  transition  region  the  occurrance  of  a 
photoelectric  edge  at  ~  11  ev  introduces  a  small  deviation  which  averages  out 
in  the  black  body  curve.  Ahead  of  the  chock,  of  course,  the  radiation  is 
highly  nonthemal  and  anisotropic  but  with  intensity  which  is  not  large  enough, 
compared  with  the  subsequent  history,  to  be  significant. 

This  conclusion  is  based  on  phenomena  encountered  at  90,000  ft  and  rather 
small  target  radii.  While  appreciably  larger  radii  are  not  of  interest  at 
this  altitude,  interest  in  larger  radii  at  higher  altitudes  suggests  that 
the  question  be  kept  open.  Indeed,  the  greater  transparency  of  the  fireball 
at  higher  altitudes  indicates  that  deviations  from  thermality  will  be  more 
important  there. 

Introduction  of  the  target  and  its  evolving  vapor  into  the  fireball  will 
perturb  to  sane  degree  the  surrounding  air  by  locally  reducing  the  tempera¬ 
ture.  The  assumption  of  thermal  radiation  at  the  unperturbed  fireball  tempera¬ 
ture,  consequently,  is  somewhat  of  an  overestimate  of  the  actual  flux  from 
the  locally  depleted  air,  Estimates  have  been  made  of  the  following  effects 
which  replenish  the  air  cooled  by  the  black  body  emission;  first,  the  air 
may  be  heated  again  by  radiative  diffusion  from  more  distant  (but  still 
nearby)  unperturbed  regions  and  second,  undepleted  air  may  be  convected  to  the 
target  by  the  streaming  velocity  induced  by  the  shock  wave. 

Assuming  that  the  local  reduction  in  temperature  is  small,  the  radiative 

diffusion  equation  may  be  linearized  so  that  the  problem  can  be  solved  in 

closed  form.  The  R.V.  is  assumed  to  be  a  perfectly  absorbing  sphere  of 

radius  a,  bathed  by  a  medium  whose  initial  temperature  is  the  unperturbed 

local  temperature  9q  of  the  fireball.  The  air  at  the  surface  of  the  R.V. 

is  assumed  to  radiate  like  a  black  body  having  the  air  surface  temperature, 

9  .  Asymptotically  at  late  time  the  solution  approaches  a  steady  state 
a 

having  dependence  on  radial  distance  r  frcm  the  R.V*  as 


1  - 


ha 


~  (ha  +  1) 


where  h  = 


and  is  the  Rosseland  mean  free  path.  At  the  surface  the  temperature 
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depends  on  time  t  as 


where 


9 

a 

9 

o 


x 


,  2 
1  "  ha" 1  ^  "  gX  erfcM]  > 

(ha  +  lj^t/a2  . 


wo  4 

The  parameter  K  =  -  contains  the  black  body  emission  flux  a9Q  and  the 

energy  per  unit  volume  E  of  the  air  at  the  unperturbed  air  temperature*  A 

time  T,  characteristic  of  the  decrease  in  surface  temperature  from  the 

initial  value  9^,to  the  asymptotic  value  Qo/(ha+l)>  is  given  by 


2 

^  _  ha  E _ 

(ha+l)  aOQ 
9 

Numerical  evaluation  of  T  and  ■■■■  shows  that  rather  different  conditions 

ha+l 

exist  in  the  shock  and  radiative  plateaus*  In  the  former,  the  local  cooling 

may  be  quite  appreciable  (for  90,000  in  the  shock  plateau  at  range 
4 

2.7  x  10  cm*  9  drops  to  .3  9  and  9  =1.5  ev)  but  t  is  rather  long 

-4  s  o  o 

(t  =  3  x  10  sec  for  the  case  above)*  Since  the  shock  region  contains 

large  air  velocities,  which  may  be  supplemented  by  the  R*V*  velocity, 

considerable  convection  occurs  during  this  time,  sufficient  to  increase 

markedly  the  asymptotic  value  (the  displacement  d  =  vt  is  d  =  3^0  cm  with 

v  -  1.1  x  10  cm/sec.  This  is  considerably  greater  than  the  assumed  R.V. 

radius  a  »  50  cm  and  much  greater  than  4_.  «  4  cm)  ,  In  the  radiative  plateau 

K 

the  mean  free  path  is  longer  principally  because  of  decreased  density  so  that 
the  asymptote  is  much  higher.  On  the  other  hand,  the  time  T  is  so  short 
that  convection  is  of  no  importance. (For  90>000  ft  in  the  radiative  plateau 
at  range  2.5  x  10^  cm,  =  *88  9q  and  9q  =  5  ev.  T  =  1*2  x  10  sec,  giving 
d  =  1  cm.  Here  a9^  v  is  reduced  to  0*6  of  the  unperturbed  medium  value*) 

For  the  cases  investigated  the  effect  of  local  cooling  is  neither  larger 
nor  completely  negligible*  Somewhat  surprisingly,  the  effect  is  smaller  at 
the  6o,OOG  ft  altitude  where  strong  hydrodynamics  decreases  the  density 
sharply  behind  the  radiative  front  to  a  value  actually  lower  than  at  a 
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comparable  stage  at  9 0,000  ft.  At  altitudes  considerably  higher  than  9°>000 
ft,  one  would  expect  that  the  decrease  of  air  density  and  increase  of  mean 
free  path  would  reduce  the  effect  to  negligible  proportions.  At  inte mediate 
altitudes  it  appears  to  be  necsssary  to  take  into  account  some  reduction  of 
flux  due  to  this  effect . 

C.  HYDRODYNAMIC  ENVIRONMENT 

As  the  shock  wave  sweeps  over  the  RoV.  the  Mach  number  of  the  air  is 
initially  large  (M1  ~  3)  and  the  stagnation  point  pressure  against  the  R.V. 
is  greater  than  ambient  by  a  larger  factor  (~  10)  .  As  the  radiation  front 
engulfs  the  RoV.  the  air  velocity  drops  and  the  sound  speed  decreases  to 
give  a  typically  subsonic  flow  (M^  ~  .7)  and  a  small  increase  of  stagnation 
pressure  over  ambient  (~  1.3)  *  The  pressure  changes  in  the  opposite  sense 
on  the  surface  away  frcm  the  explosion  center,  of  course.  The  velocity  of 
the  R.  V.  itself  is  not  negligible  compared  with  the  air  velocity  and  may 
add  to  or  subtract  from  that  of  the  air  depending  on  the  relative  orientation 
of  the  R.V.  in  the  explosion.  While  these  dynamic  effects  are  not  small,  the 
sign  of  the  effect  depends  on  R.V.  orientation  and  surface  of  interest.  The 
ambient  pressure  is,  in  a  sense,  the  central  value  of  the  pressure. 
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III.  NATURE  OF  THE  IMEHACTION 

Calculations  carried  out  with  the  code  as  described  in  Section  I  and 
with  boundary  conditions  for  the  three  cases  recorded  in  Section  II. A.  are 
analyzed  in  this  section.  The  properties  assumed  for  the  surface  of  the  R.V. 
are  those  characteristic:  of  a  low  thermal  conductivity  plastic  composed  of 
low-2  elements  for  which  more  data  are  recorded  in  Section  I.G.  For  the 
purpose  of  taking  into  account  the  divergence  of  the  expanding  vapor  as 
described  in  Section  I.C.  a  characteristic  R.V.  radius  r  =  50  cm  was  used. 

In  each  case  an  outer  layer  of  the  plastic  was  divided  into  150  pson^s,'  the 
outermost  one  of  which  was  10  gm/cm  in  thickness.  Zones  increased  in 
thickness  with  increasing  depth  into  the  surface  to  give  a  total  zoned  layer 
thickness  of  0.5  gm/cnT,  or  C?5  cm  at  the  initial  density  of  1.0  gm/cm^. 

-3 

Calculations  were  continued  forward  in  time  until  from  1»5  x  10  to 
-3 

4  x  10  sec  had  elapsed  -  This  history  corresponds  to  transit  through  the 
shock  plateau  and  rise  of  source  temperature  toward  its  peak,  stopping 
somewhat  before  achieving  the  peak  source  temperature.  Impulse  experienced 

by  the  R.V.  in  each  of  these  cases  by  the  time  of  problem  termination  was 

5  2 

approximately  lCr  dyne  -  sec /cni  . 

Each  problem  required  typically  a  little  more  than  2000  cycles  or  time 
steps  and  a  little  less  than  one  hour  of  7090  computer  use.  Calculations 
were  re  zoned  to  remove  superfluous  zones  as  often  as  twice  per  problem  and 
in  one  of  the  problems  (h  =  90,000  ft,  r  =  2-5  x  10^  cm)  almost  a31  of  the 
mass  provided  in  the  mesh  had  been  ablated  by  time  of  termination. 

Results  of  the  calculations  are  examined  in  the  sections  below.  In 

Section  A  for  one  of  the  problems  the  spatial  characteristics  of  temperature 

0,  pressure  P,  velocity  v,  density  p,  and  Rosseland  mean  free  path  are 

ft 

examined.  A  comparison  between  the  radiation  fluxes  calculated  with  and 
without  frequency  dependence  Is  also  made.  In  Section  B  all  three  problems 
are  analyzed  for  time  dependence  of  impulse  per  unit  area,  l/A,  and  ablated 
mass  per  unit  area,  m&.  The  coupling  constant,  t/§,  is  also  examined  during 
the  interaction  history.  Section  C  contains  a  comparison  of  two  problems 


SFHRFT 


SECRET 


26 

differing  only  in  the  way  source  energy  is  deposited  in  vapor  and  solids . 
Finally,  in  Section  D  the  conditions  for  convective  instability  are  examined 
and  estimates  are  made  of  the  role  it  plays  in  the  TME. 

A.  SPATIAL  CHARACTERISTICS 

In  Figs.  III-l  to  III -20  profiles  of  five  quantities  are  shown  at  four 

4 

instants  during  the  90,000-ft  explosion  at  interaction  range  =  2.5  x  10  cm. 
In  Fig.  III-l  the  temperature  at  t  =  4.3  x  10 sec  is  seen  to  be  almost 
constant  at  0.45  ev.  At  this  time  energy  from  the  fireball  streams  through 
the  cool  vapor  with  very  little  absorption  so  that  most  of  the  source  is 
deposited  in  the  solid  where  it  is  very  efficient  in  producing  vapor.  The 
boundary  pressure  at  the  interface  is  important  in  retarding  the  blowoff  and 
in  elevating  the  pressure  on  the  solid  above  the  free  blowoff  value.  Also 
for  this  reason  the  blowoff  velocity  remains  low  and  the  vapor  density 
relatively  high.  The  limitation  of  the  temperature  at  0.45  ev  is  due  to  a 
combination  of  the  large  dissociation  energy  needed  to  reach  higher  tempera¬ 
tures  and  the  small  absorptivity  of  the  vapor  at  low  temperature.  The  energy 

r  -4 

absorbed  from  the  source  does,  however,  at  t  ~  6.3  x  10  sec  exceed  the 
dissociation  energy  in  the  outside  layers  resulting  in  a  marked  temperature 
rise  as  shown  in  Fig.  111-6.  Accompanying  the  increase  in  temperature  is  an 
increase  in  vapor  ionization  and  absorptivity  which  markedly  affects  the 
behavior  of  the  problem.  The  solid  suddenly  is  shielded  from  the  source  by 
the  very  opaque  vapor  which,  in  turn,  is  heated  to  a  high  enough  temperature 
for  reradiation  to  play  a  role.  In  Fig.  III-ll  at  8.0  x  10  ^  sec  the  vapor 
temperature  has  reached  3-9  ev,  sufficient  to  reradiate  35$  of  the  incident 
energy  from  the  system.  The  source  is  now  being  absorbed  by  the  denser  vapor 
at  the  head  of  the  advancing  temperature  wave  which  is  also  radiating  energy 
ahead  into  the  cold  vapor.  Considerable  hydrodynamic  disturbance  results 
from  the  rather  sudden  expansion  of  the  previously  cold  vapor  and  impels  a 
compression  and  shock  wave  forward.  The  hot  vapor  now  expands  outward  with 
velocity  ~  10  but  the  cold  vapor  is  recompressed  against  the  solid. 

In  Fig.  Ill- 16  the  temperature  wave  has  progressed  inward  close  to  an 
equilibrium  position  in  the  vapor.  Reradiation  does  not  yet  contribute 
appreciably  to  the  ablation  rate  but  will  do  so  at  a  little  later  time.  It 
is  expected  that  on  completion  of  the  vapor  heating  a  greater  fraction  of 
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source  energy  will  reach  the  solid  which,  together  with  reradiation,  will 

appreciably  increase  the  ablation  rate.  Loss  of  energy  fran  the  system  by 

reradiation  still  accounts  for  38$  of  the  source.  At  this  time  the  inward  - 

7 

driven  wave  has  reflected  from  the  R.V.  surface  with  pressure  p  =  9  x  10 
and  is  driving  the  vapor  outward  with  v  >  2  x  10  cm/sec. 

Throughout  these  problems  the  boundary  pressure  has  ensured  that  the 
density  remain  high  enough  for  the  hot  vapor  to  stay  in  the  radiation  diffusion 
regime  o  As  a  consequence  of  this,  local  thermodyiiaraic  equilibrium  has  been 
maintained  throughout  the  entire  calculation. 

One  calculation  of  the  frequency-dependent  radiation  transport  was  also 
made  on  the  above  problem.  Corresponding  to  t  =  8.5  x  10  sec,  the  calcula¬ 
tion  permits  comparison  with  the  transport  calculation  using  the  Rosseland 
mean  free  path  as  employed  for  the  entire  histories  of  the  time -dependent 
calculations  *  At  this  time,  as  noted  above,  the  temperature  profile  consists 
of  a  high  temperature  plateau,  an  abrupt  transition  region,  and  a  low- 
temperature  plateau  into  which  radiation  streams  from  the  hot  region.  In 
Pig.  III-21  the  radiative  fluxes  calculated  in  these  two  ways  are  shown 
throughout  the  system.  At  the  outer  boundary  the  flux  is  positive,  corres¬ 
ponding  to  black  body  emission  frem  the  high  temperature  plateau.  Deeper 
into  the  vapor  the  flux  decreases  and  changes  sign  as  the  region  of  transi¬ 
tion  is  approached.  This  change  of  sign  of  flux  is  noted  on  the  graph  by  a 
vertical  line  marking  the  position.  In  the  negative  flux  region  the  flux 
reaches  a  minimum  as  the  temperature  falls  and  then  rises  again  toward  zero 
by  absorption  toward  the  solid.  Good  agreement  is  found  in  the  high  tempera¬ 
ture  regime  •  In  the  transition  region,  however,  the  Rosseland  curve  is  too 
high  while  at  the  solid  it  is  appreciably  too  low.  In  view  of  the  smallness 
of  the  radiative  flux  compared  to  the  source  flux,  very  little  error  was 
made,  showing  that  the  Rosseland  treatment  is,  in  this  case,  quite  satisfactory. 
For  times  when  radiation  contributes  importantly  to  the  ablation  rate,  more 
important  deviation  may  be  expected,  however.  Work  now  in  progress  on  using 
the  frequency-dependent  code  throughout  the  calculation  may  then  be  very 
valuable.  In  Fig.  III-22  the  frequency  dependence  of  the  flux  is  displayed 
at  several  positions  in  the  vapor. 
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B.  EVOLUTION  IN  TIME,  LETHALITY,  AND  MATERIAL  LOSS 

Luring  the  interaction  the  pressure  at  the  solid  surface  results  in  a 
continued  increase  of  impulse  contributed  to  by  both  the  ambient  pressure 
and  the  ablation  pressure  as  influenced  by  the  hydrodynamics  of  the  vapor. 

The  amount  of  mass  ablated  determined  by  source  and  re  radiation  also  increases 
with  time.  In  Figs.  III-21  through  III-23  the  history  of  both  of  these 
quantities  is  shown  for  the  three  interaction  problems.  At  early  times  in 
the  shock  plateau  the  curves  are  all-  very  similar  with  approximately  linearly 
increasing  ablation  and  impulse  with  time.  This  phase  is  terminated  by  the 
progressive  increase  in  source  rate  as  the  radiative  wave  arrives.  Happening 
at  earliest  time  in  Problem  1  and  successively  in  Problems  2  and  3,  the 
ablation  rate  responds  immediately,  followed  much  more  slowly  by  the  impulse 
curve.  Problem  3  was  terminated  in  this  state.  The  first  two  of  the  problems* 
however,  show  the  effects  of  the  catastrophic  surface  temperature  increase 
reflected  by  sudden  decrease  in  the  mass  ablation  rate  and  in  the  impulse 
curve  by  a  sudden  drop  in  the  ablation  pressure,  introducing  a  plateau  in 
the  impulse  curve. 

In  Problem  2  the  calculation  was  ended  almost  immediately  after  ioniza¬ 
tion  started.  Only  in  Problem  1  is  this  phenomenon  followed  farther.  The 
mass  ablation  curve  is  greatly  decreased  in  slope  for  an  extended  interval 
equal  to  at  least  10  sec  for  Problem  1.  In  the  impulse  curve,  however, 
structure  is  introduced  by  the  violent  hydrodynamic  motions  triggered  by  the 
reduction  in  ablation.  Pressure  at  the  solid  interface  is  reduced  at  the 
same  time  that  a  shock  wave  is  driven  inward  by  the  expanding  material  in 
the  heating  wave.  These  then  cooperate  to  raise  the  pressure  of  the  solid 

Q 

momentarily  to  P  >  10°  as  the  shock  reflects.  Consequently,  the  impulse 
curve  rises  abruptly  and,  after  some  smaller  fluctuations,  assumes  a  new 
slope  smaller  than  the  original.  While  Problem  1  was  terminated  before  the 
thermal  wave  reached  the  solid,  this  evert  should  initiate  a  new  increase  in 
both  mass  ablation  rate  and  ablation  pressure. 

It  is  also  of  interest  to  examine  the  efficiency  of  the  thermal  energy 
source  in  producing  impulse.  This  has  been  done  by  forming  I/E  =  d(l//S)/odt 
where  Q  is  the  thermal  power  per  unit  area  from  the  source  Q  -  in 

watts/cm2.  In  Figs.  111-24  through  111-26,  l/E,  the  instantaneous  coupling 
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coefficient,  is  shown  vs  time  for  the  three  problems  =  The  values  are  sur¬ 
prisingly  large  at  early  times  in  consequence  of  the  fireball  ambient 
pressure.  This  pressure,  more  conventionally  thought  of  as  part  of  the  blast 
wave  and  associated  with  the  blast  regime,  is  relatively  more  important  in 
the  shock  plateau  before  the  ablation  rate  has  risen .  At  the  60,000-ft 
altitude  the  interaction  is  almost  entirely  attributable  to  the  blast  regime 
with  only  minor  TME  superimposed.  This  altitude  is  clearly  below  separation 
boundary  between  blast  and  THE.  At  the  altitude  of  90>000  ft,  however,  blast 
and  thermal  effects  are  interwoven  into  a  closely  interrelated  cumulative 
effect.  The  more  specifically  ablative  values  of  i/E  are  obtained  as  the 
ablation  pressure  rises  well  above  ambient .  While  considerable  structure  is 
introduced  into  the  Problem  1  curve  by  vapor  heating,  typical  values  of  i/E 
associated  with  Lhe  thermal  mechanical  effect  appear  to  be  1-5* 

Lethality.  It  is  striking  that,  as  a  consequence  of  vapor  shielding 
Problem  1  at  a  shorter  range  from  the  explosion  has  less  ablation  than  and 
comparable  impulse  to  Problem  2  after  1.5  x  10  ^  sec  of  interaction.  The 
impulse  curve  for  Problem  2  is  considerably  steeper  and  thus  more  damaging, 

«  Ithough  both  locations  have  probably  suffered  lethal  damage.  (Lethality 

is  evaluated  in  terms  of  criteria  established  by  Abrahamson  and  Lindberg^^ ^ 

4  ,  7  x 

These  require  total  impulse  per  unit  area  >  2  x  10  and  pressure  >  4  x  10  „ ) 

The  vapor  shock  in  Problem  1,  however,  is  worthy  of  further  investigation 

as  to  its  lethality.  The  impulse  curve  for  h  =  60,000  ft  in  Problem  3  is 

considerably  shallower  corresponding  to  lower  pressure  which  is  probably 

not  lethal  to  a  hardened  R.V. 

Material  Loss.  The  amount  of  material  lost  by  ablation  is  not  large 
during  the  calculations.  If  the  rate  of  ablation  were  maintained  near  its 
peak  value  for  the  duration  of  the  high  temperature,  however,  significant 
material  loss  would  result. 

C.  COMPARISON  OF  SOURCE  FUNCTIONS 

While  for  an  isotropic  incident  distribution  the  correct  law  of  decrease 
of  flux  at  optical  depth  x  below  the  surface  is  given  by  PE^fx)  (see  Section 
I. A.  for  a  description  of  the  source  subroutine),  an  early  versicn  of  the 
source  erroneously  contained  E?(x).  While  these  functions  lead  to  the  same 
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total  incident  flux,  the  latter  function  deposits  the  energy  nearer  the 

k 

surface  than  the  former.  Since  some  cf  problem  h  =  90,000,  r  =  2.5  x  10 
was  run  with  E^  before  being  repeated  with  a  comparison  exists  between 

problems  differing  only,  in  the  angular  distribution  of  the  source  or 
equivalently  in  the  deposition  function  for  the  source  energy.  This  compari¬ 
son  gives  an  indication  of  the  sensitivity  to  depth  of  deposition  of  results 
and,  consequently,  gives  some  basis  for  judgement  of  how  accurately  the  non¬ 
isotropic  distribution  at  higher  altitudes  must  be  approximated.  In  Fig. 

Ill -27  the  impulse  per  unit  area  and  mass  ablated  are  compared  at  rather 
early  times.  While  initially  the  curves  are  substantially  the  same  with  all 
of  the  energy  penetrating  to  the  solid,  as  vapor  is  evolved  the  function 
deposits  more  energy  in  vapor,  leaving  less  for  vaporization.  The  greater 
deposition  in  the  vapor  at  the  air  interface  supplies  the  dissociation  energy 
there  at  an  appreciably  earlier  time  and  initiates  the  catastrophic  heating 
and  ionic  absorption  which,  in  turn,  is  the  reason  for  the  reduction  in  rate 
of  increase  of  impulse  and  mass  ablation.  Since  these  phenomena  occur  at 
earlier  time,  source  rate,  and  total  energy  the  deposition  reduces  the 
interaction  considerably.  It  is  not  known,  of  course,  whether  the  earlier 
completion  of  the  vapor  heating  would  result  in  a  compensating  earlier 
resumption  of  strong  ablation.  The  differences  in  Fig.  III-27,  however, 
occur  at  a  significant  part  of  the  interaction  in  that  the  impulse  is  increas- 
ing  most  rapidly  there  and  are  of  sufficient  magnitude  to  indicate  considerable 
sensitivity  of  results  to  angular  distribution.  Some  care  in  specifying  the 
distribution  is  indicated.  The  same  may  be  said  of  the  low  temperature  vapor 
absorption  coefficient,  pointing  up  the  need  for  more  accurately  known 
material  properties. 

D.  CONVECTIVE  INSTABILITY 

There  is  theoretical  evidence  that  radiation  is  not  the  only  energy 
transport  mechanism  present  in  the  vapor  when  the  temperature  reaches  values 
large  enough  to  make  radiation  itself  an  Important  process.  The  additional 
transport  mechanism  is  one  which  takes  place  if  elements  (bubbles)  in  the 
gas  may  be  convected  by  the  acceleration  field  to  increasingly  greater 
distances  from  their  point  of  origin.  Only  if  the  motion  is  unstable  does 
the  convection  establish  a  non-zero  steady  state  in  which  energy  is  convected. 


QFPPFT 

\J  I-  \J  l  \  L. 
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.An  analog  of  this  process  is  convective  heat  transfer  in  the  atmosphere  where 
the  buoyancy  displacement  of  air  cells  may  under  certain  conditions  be 
unstable.  There  is  experimental  evidence  from  the  Orion  investigations 
indirectly  applicable  to  the  thermal  interaction,  that  the  process  is  present 
in  a  similar  radiative  ablation  environment. 

Consider  a  slab  of  material  with  a  positive  pressure  gradient  and  with 
a  density  and  temperature  gradient  as  shown  below. 


P 


The  two  curves  labeled  0  and  p  indicate  the  way  in  which  the  tempera- 

a  a 

ture  and  density  would  vary  if  from  point  A  the  flew  were  adiabatic  and  the 
pressure  remained  the  same.  If  a  perfect  gas  equation  of  state  is  assumed, 

P  =  op9,  then  along  P  =  Kp^,  the  adiabat  is  determined  by  P: 


P 


a 


(?)1/y 

V 


M 

dzl  a 


£_  dP 
YP  dz 


and 


9 

a 


k  ,p1“  y  •  a©  j 
5  '■K;  *  "  dzl  a 


9  dP 
P  dz  ’ 


A  small  element  of  vapor  displaced  to  the  left  from  A  along  the  adiabat 
experiences  an  unbalanced  force  on  the  element  proportional  to 


If  the  quantity  given  in  the  brackets  is  negative,  then  the  element  will 
continue  to  move  to  the  left  and  the  vapor  will  be  unstable  to  small  pertur¬ 
bations.  An  element  displaced  to  the  right  from  B  will  be  unstable  with 
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respect  to  displacement  to  the  right  under  the  same  conditions.  The  result 
is  convective  energy  transfer  toward  the  solid  vapor  interface  due  to  the 
circulation  in  the  unstable  layer  of  hot  elements  toward  the  colder  region 
and  cold  elements  toward  the  hotter  region.  The  condition  for  instability 
can  be  written  as 


BP  /» 

Bz  '  Bz  a 


0  . 


Convective  energy  transport  in  an  unstable  layer  in  the  presence  of 

(12) 

energy  transport  by  radiation  is  also  well  known  in  stellar  interiors. 

A  similar  energy  flow  has  been  observed  in  fluids  in  the  presence  of  thermal 
(13) 

heat  conduction'  where  the  convective  mixing  cells  that  are  formed  are 
known  as  Benard  cells. 


The  energy  flux  in  the  unstable  layer  can  be  estimated  in  two  ways. 


(12) 

Schwarzs  child'  'derives  a  convection  flux  that  is  given  by 


</>  = 
rc 


PC 


■  L- 


BP 

pBz 


r 


where 


.  / _  BO  B9 
^z}  ~  Bz  a  "  Bz 


and  l  is  the  mixing  length.  This  quantity  can  be  estimated  to  be  the  width 
of  the  unstable  layer. 

In  Kef.  13  the  effect  of  viscosity  in  the  Benard  cell  problem  is 
considered.  It  is  shown  that  the  energy  flux  is  unchanged  until  a  critical 
value  for  a  dimensionless  number,  the  Rayleigh  number  (R),  is  reached  in  the 
layer.  A  Rayleigh  number  can  be  defined  for  the  radiation  environment  as 


R  = 


.  9P  .  Sfi.  .  1  U)b 
Bz  Bz  Bz  T| 


where  is  the  radiative  flux  and  Tj  is  the  viscosity.  The  ratio  of  energy 

transfer  with  the  instability  present  to  the  energy  transfer  in  the  absence 
of  convective  mixing,  the  Nusselt  number  (Nu),  has  been  measured  as  a  function 
of  R  for  the  Benard  experiment.  The  results  are  independent  of  the  fluid  used 
in  the  experiment.  The  critical  value  of  R  for  the  experiment  reported  is 
~  l60O.  The  Nusselt  number  is  observed  to  be  a  very  weak  function  of  R, 
e.g.,  for  R  =*  10^,  Nu  --  8. 
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Fig.  III-10 
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Fig.  111-15 
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APPENDIX  A 

THERMAL  SOURCE  FLOW  DIAGRAM  AND  CODE 

The  quantities  required  as  input  for  the  "QUE"  and  ,lHYDRO"  subroutine 
are  given  below,  together  with  descriptions. 

HYDRO 

P(-j)  Fractional  weights  of  a  finite  number  of  frequency  groups  for  the 

source . 

CRTPC(j)  Average  energies  of  a  finite  number  of  frequency  groups  for  the 
source. 

CHRN0(24)  Number  of  frequency  groups  of  the  source. 

CHRN0(25)  Constant  coefficient  used  for  the  CAPAC(l)  calculation. 

RDK(l-52)  Times  (t.)  for  the  QUE  subroutine. 

X 

RDK( 53-104)  Initial  boundary  temperatures  (9  )  for  the  QUE  subroutine. 

Li 

SOLID(l)  Positive,  zero  or  negative  trigger  to  use  EXPi4,  EXPI3;  or  EXPF 
in  source  angular  distribution. 

S0LID(2)  Cosine  of  angle  of  incidence  of  source  ray,  usually  one. 
S0LID(3)  Constant  factor,  £  1,  to  represent  a  distant  thermal  source. 
S0LID(8)  Initial  source  rate  used  to  establish  conduction  wave. 

S0LID(9)  Constant  coefficient  used  for  the  XKLC  calculation. 

A  Atomic  number  of  target  material. 

Z  Charge  number  of  target  material. 

QUE 

THETAK(l-52)  Times  (t.)  for  the  pressures  in  "hydro". 

THETAK( 53-104)  Boundary  pressures  (9^  in  "hydro". 
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pOFR(lB)= 


Normalized  l'lux  at  surface  Is  1. 


IK 


-  C1CIT.‘0(?H )  |  Load  the  number  of  frequency  groups. 


pb(j)  »  •;  | 


initially  5.:et  each  frequency  group 
optical  depth  to  zero. 


|  I  =  IBM1  | 

r  Boih 

h 

< 

> 

L _ _ 

Test  whether  a  zone  is  boiling. 


TEMP(h)  -  3V(  IAM  1) 
SV(lhMl)  *  SV(l/iMl  -  1) 
NBOTT.  =  1 

A( X )  =  A(I)/(1.-S0LID(2^) 


Boiling  zone  source  is  being  calculated  - 
Get  up  quantities  to  evaluate  optical 
depth  of  unvaporized  fraction  of  the 
boiling  zone . 

Exit  to  code  to  form  contribution  to 
optical  depth  at 
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DO  -x 

CJ=1,  IK 


CRTPC(j)o 

_ O 

_eii) _ 


u  - 

uJ: 

TH-1P(S)=  (eu.i)+  o.2 


TEMP(2)  =  ru)eUk+0.2 


Form  frequency-dependent  part) 
of  absorption  coefficient  and 
cum  all  f'roup  contributions  ! 
to  optical  depth  and  flux.  1 
Linear  interpolation  between  ' 
low  temperature  and  ionic  * 
coefficients  when  <Z>  <  .>  1 


----- 

> 

[rE3lP(F)  -  (<Z>  •  +  (.></.»•  XKLO)  •  2.0  ~~] 


PB(J)  -  v  pb(j) 


TEMP(  2  )=S0LID(  2  )* 
EXPF(-PB(J)/S0LID(2)) 


^  GO  TO  410 


T1KP(  2 )  =EXPI 3  ( PB(  I ) ) 
TKMP(  2 )  =T>JMP(  2 ) +TKM  P(  2  ) 

GO  TO  410 


T1MP(2)=3.*EXPI4(PB(J))  | 


|G0FR(l)  -  l?34P(2)*D(j)+  GOrH(l) 
410 
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Determine  whether  a  boiling  zone 
is  being  calculated  and  return 
to  proper  location. 


0 

_ v _ _ _ 

[lO'fo  (jjj.iQ,  .'ROIL 
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I L=SUMRR 
DO  1  I  =  I L  ,  IBM1 
1  GOFR ( I )=0o 
SUMRB=  I R 

C  T{I)  GREATER  THAN  OR  EQUAL  TO  ZERO 

C  T1I  +  1!  GREATER  THAN  OR  EQUAL  TO  T 1 1  ) 

IF  ( RDK ( 1 ) )  180,110,110 
110  IF  (RDK(l)-TH)  120,120,161 
120  DO  135  J=2,52 
J  =  J 

IF  (ROM  Jl-RDKl  J-l))  170,170,130 
130  IF  (RDK(J)-TH)  135,150,140 
135  CONTINUL 
GO  TO  160 

140  TEMP ll)  =  RDK ( J+51)+(RDK( J  +  52 ) -RDK ( J+ 5 1 )  ) / (ROK < J J-RDK ( J-l )  )* 
1 (TH-KDK( J-l) ) 

GO  TO  200 

150  TEMP(1)=R0K( J+52) 

GO  TO  200 

160  S6=0„ 

161  S0Lin(8)=0o 

GO  TO  1000 

170  IF  ( RDK ( J ) )  180,160,180 

180  WRITE  OUTPUT  TAPE  6 , 190 , RDK { J ) , J 

190  FOR  MAT ( 1 9H1  INPUT  THETA  ERROR , 1 P El  5 . 7 , I  4  ) 

CALL  UNCLE 

200  50Lini3)  =  l„026E12*SULIDI 3)*TEMP<  1!**4 
GOFR ( I B ) = 1 o 
I K=CHRNO ( 24 ) 

DO  30  J=1 ,IK 

30  P  b  !  J  )  -=  0  o 
I  =  I  B  M 1 

4  I F ( BQI LA ) 31 , 3 , 31 

31  IF ( I -I  AMI )3»  32 ,3 

32  TFMP(4)=sSVUAMl) 

SVt I  AMI ) =  SV(  I  AMI-1 ) 

NBD I L= 1 

All )=A(I) / (1 O-S0LID124)  ) 

GO  ro  3 

33  SVt IAM1 )  -  SOL 10(27) 

NBOI L=2 

All  )=AU  } /SOLID! 24 )*(1«-S0LID124)  ) 

GOFR  (  I  '  =0  . 

GO  TO  3 

37  A { I ) -A ( I } *SOLID( 24) 

S  V i I  AMI ) =TEM P { 4 ) 

5  SMLQl  I  }  ■=  l  GOFR l  I  +  1)-G0FR(  I)  )*S0LID(8) 

IFIGOFRI I >-.0001)10, 27, 27 

27  IFl [-1)29,29,9 

29  SOLID ( 8 ) =SGL ID ( 8 ) * ( 1 •— GOFR ( I)) 
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GO  TO  28 
9  1=1-1 
GO  TO  4 

10  SMLCM I )=SMLQ( 1 ) +GOFR ( I ) *SOL ID( 8) 

28  SUMR8-- 1 

RPltt=RPI B+SOLIOl 8 ) *DTH2 
1000  RE TURN 

3  XKLC  =  SOL  10(9 )  /SV ( I ) 

I  F  (  "  £ W ( I  )  )  511,511, 509 

511  TEMP ; 11 ) =XKLC/AI I ) *G  !  I  ) 

DO  512  J= 1 , 1  K 

PB(  J)=PR(  J)+TEMP(ll  ) 

IF  (  SOL  Il)(  1) )  420,421  ,422 

420  T  E  M  P  '  2 ' -  SOLI  0(25  »EXPF ( -PB( J ) /SOLID! 2 ) ) 

GO  TO  '>]  2 

421  TEMP ! 2 !=£XPI3(PR!J)  ) 

TEMPI/.  I  =  TEMP ( 2 ) +  TEMP ( 2  ) 

GO  1!)  512 

422  TEMP (2 5=3. *EXPI4(PB( J)  i 

512  GOFR ; I )= TFMP ( 2 )»D< J ) +GOFR( I ) 

GO  TO  513 

509  00  610  K=l, MAXIM 

I F ( l-LMOAIK+l 5)510,510,610 

510  1 F  =  K 

GO  TO  630 
610  CONTINUE 

63  0  TLMP(  135  =  13.  6*CHRN0  (  IF  }**?/  THETAt  I  ) 

TEMP ( 14 ) = THE  TA(I)**1.5*SV( I > /FEW ( I) *1 . E 10/ AMASNO ( I F ! 
TEMP ( 1 5 ) =LOGF ( TEMP ( 14) ) 

TEMP (12)=MIM1F (TEMP ( 13 ) , TEMP (15 ) ) 

C  APAC  ( I  )  -  CFIR  Nu  i  2  5  )  *  FEW  (  1)  »M  AX1F 1 1 .  ,  FEW  (  I  )) 

1*AMASN0( IF)/ (THETAt I) * *2* TEMPI  14) J 
DO  40  J  =  1 » J K 

T£  MP  ( 1 1 )  ---CRT  PC  (  J  )  *TE  MP  ( 1 )/ THETA  (  I  ) 

IF (TEMP (13 ) -TEMP (12) 5800,950,950 
800  TCMP(2 )-CAPAC( I )/TEMP(ll )**3*(EXPF(TEMP( 11) )-l.)+„2 
GO  TO  35 

95  0  TEMP ( 2 ) =C APAC ! I  5 /TEMPI  11 )**3*( l.-EXPF(-T  EMP  (  11 5  5  )* 

IE  XPF  {  TEMP  (1?  )  H-.2 
35  I  FMP  {  I  i  ;  -  „5  FEW  (  I) 

IF! TEMP f 11) ) 39,39,38 

38  TEMPI  2'.  =  (  FEW  C I  )  *TEMP  (  2  ) +T  EMP  ( 1 1  )*XKLC)*2  . 

3  9  PB(J)=TEMP(2)/A(H*G(I  5  +PR  !  J  ) 

IF(SOLl(!(l)!  400,401,402 

400  TEMP (2 i=S0EIU(2)*EXPF(-PB( J ) /SOL  ID (  2 ) ) 

GO  TO  410 

401  T C M P  C  2 ) =  E XPT  3 ! P 8  ( J  )  > 

TEMP ( ? ) - T EMP ;  ? 1 +TEMP ( 2  ) 

GU  n:  410 

402  TEMP  (2) -*3. »L XT'  1 4 ( PG { J )  ) 
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410  G0FR{ I l-TTMP (2 ) *D ( J ) +GOFR { I ) 
40  CONTINUE 
513  IF(ROILA)  36,5,36 
36  I F ( I— I  AM 1)5,11,5 
11  GO  TO  ( 33 ,37 ) ,NHUIL 
F.  NO 
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PRESSURE  BOUNDARY  CONDITION  FLOW  DIAGRAM  AND  CODE 
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I  corr-iiai-. 
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! ((-iw w( j J-thi*1  My--)',')  - ij’oaPiO  ni!®u;(j  .3?)  j — ^Vj) 


1.6^)— ((rHf,T.\i.(j))  J 


TL7‘P(  )  :ll!.V(\r.(  J  •  Jl )  •  (THi7PAK(.;  •  '):■■)  -THi?IVK(j 1  31)) 

/(  J’ilKV  •*(.I)-,.UK*a(c -!))*(  rH-THijT.1!C{  J  - ! ) ) 
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SUBROUTINE  HYDRO 


(w) 
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IF  I  THETAM 1  )  )  180,100,100 
100  IF  ( THETAK ( 1 J-TH )  115,115,160 
115  DO  135  J=2,52 
J=J 

IF  ( THETAM  J)-THETAK(J-l) )  170,170,130 
130  IF  (THETAM J )— TH )  135,150,140 
135  CONTINUE 
GO  TO  160 

140  TEMP (l)=THETAK(J+51)+(THETAK{J+52) -THETAM  J  +  51 ))/ 

1 ( THETAK ( J)-THETAK(J-l) ) » ( TH-THETAK ( J-l ) ) 

GO  TO  10 

150  TEMP ( 1 ) =THET AK ( J+52 ) 

GO  TO  10 
160  TEMP ( 1 ) =0 . 

GO  TO  10 

170  IF  (THETAKIJ})  180,160,180 

180  WRITE  OUTPUT  TAPE  6 , 190 , THETAK ( J ), J 

190  FORMAT l 22H1  INPUT  PRESSURE  ERROR, 1PE15 . 7 , 14) 

CALL  UNCLE 

10  IF  (PUSHA)  11,20,11 

11  P ( I B ) =  TEMP ( 1 ) 

IC= I API 

C ( I  A) *R ( I  A) 

C ( I  AMI ) =R( IAM1 ) 

IF  (lBPl-IG)  12,12,14 

12  DO  13  I=I8P1,IG 

13  R(I)*R(I) +R0 ( I ) *DTH2 

14  IF  ( I B— I A )  120,120,50 
20  IC=IA 

IF  (PUSHB)  40,30,40 
40  P( IAM1 )=TEMP (1 ) 

RD( IB) =ASMLR 
RIB*Rt  XR) 

GB=G( IB) 

G(IB)=2.*(GL-S0LID(30) ) 

30  IF  ( THETAM2  ) )  180,50,31 

31  P ( I R ) =TEMP( 1 ) 

50  CONTINUE 

80  DO  110  I=IC, 18 
Cl 1 1 -R ( I ) 

RO 1 1 )  =Rt)  ( I )+  (2.*  (P(  I-l)-P(I)  )/(  G(I)+G(I-m»A(  I)  )  *DTH1 
R (I) =R (I) +RD ( I ) *UTH2 
110  CONTINUE 
120  RETURN 
END 
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APPENDIX  C 
PROPERTIES  OF  CEL 

The  DIAPHANOUS  data  are  arranged  in  ascending  values  of  the  temperature 
as  the  first  independent  variable  for  each  of  which  the  second  independent 
variable  is  the  density  arranged  in  descending  sequence*  Each  page  contains 
values  for  a  single  temperature  (ev)  and  lists  for  reference  the  integrated 
black  body  thermal  power  (ergs /cm2 sec) ,  In  the  body  of  the  table  the 

following  quantities  are  tabulated: 


Quantity 

Units 

Explanation  Formula 

RKO 

Material  density 

P 

log  JRHO 

- 

Natural  logarithm  of  density 

- 

Pressure 

ergs /cm^ 

Material  pressure 

P 

ENERGY 

ergs/gm 

Total  specific  material  energy 

E 

EION 

ergs/gra 

Ionic  specific  material  energy 

- 

ZBAR 

- 

Number  of  free  electrons  per 

atom 

EGAM 

- 

Strong  shock  "gamma"  y  = 

!  +  ~ 
pE 

KROS 

cm2 /gm 

Rosseland  mean  opacity 

kh 

KFLK 

cm2/gm 

Emission  or  Planck  mean  opacity 

HP 

M3AR0S 

cm 

Rosseland  mean  free  path  X^ 

_ 

V 

LMDAPLK 

cm 

Planck  mean  free  path  X^ 

_ 1_ 

=  V 

TCOOL 

sec 

Radiative  emission  cooling  t  = 

E 

- II — 

time 

kaQ  H 

p 
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DIAPHANOUS  DATA  CH2 

THETA ( EV )  ^  1.500  SIGMA  THETA  4  = 


GENERAL  ATOMIC 
5.194+12 

LOG  RHO  PRESSURE  ENERGY  EION  ZBAR 

(GM/CM3J  ( ERG/CM3 )  CERGS/GM)  (ERGS/GM) 


RHO 

{ GM/CM3 } 

2.340-01 

2.211-02 

1.620-03 

1.226-04 

1.202-05 

1.638- 06 

2.639- 07 
4.324-08 
6.334-09 
7.229-10 


RhU 

( GM/ CM3 ) 

2.340-01 

2.211-02 

1.620-03 

1.226-04 

1.202-05 

1.638- 06 

2.639- 07 
4.324-08 
6.334-09 
7.229-10 


-1.452+00 

-3.812+00 

-6.425+00 

-9.007+00 

-1.133+01 

-1,332+01 

-1.515+01 

-1.696+01 

-1.888+01 

-2.105+01 


LOG  RHO 
(GM/CM3) 

-1.452+00 
-3.812+00 
-6.425+00 
-9.007+00 
-1.133+01 
-1.332+01 
-1.515+01 
-1.696+01 
-1.888+01 
-2. 105+01 


7.336+10 

7.073+09 

5.521+08 

4.842+07 

5.837+06 

9.315+05 

1.598+05 

2.667+04 

4.005+03 

4.8^7+02 


KROS 

CCM2/GM) 

1.036+04 
9.107+03 
8.215+03 
5.551+03 
2.448+03 
6.646+02 
1.308+02 
2.545+01 
6.575+00 
i. 738+00 


5.748+11 

6.352+11 

8.235+11 

1.320+12 

2.207+12 

3.063+12 

3.453+12 

3.593+12 

3.861+12 

4.669+12 


KPLK 

(CM2/GM) 

3.159+05 

3.056+05 

2.868+05 

2.356+05 

1.429+05 

5.278+04 

1.271+04 

3.048+03 

1.137+03 

4.284+02 


1.046+11 

1.553+11 

3.123+11 

7.278+11 

1.479+12 

2.210+12 

2.555+12 

2.668+12 

2.913+12 

3.653+12 


LMDAROS 

(CM) 

4.125-04 

4.966-03 

7.514-02 

1.469+00 

3.398+01 

9.188+02 

2.896+04 

9.088+05 

2.401+07 

7.960+08 


0.019 

0.039 

0.107 

0.283 

0.578 

0.848 

0.967 

1.004 

1.054 

1.201 


LMDAPLK 

(CM) 

1.353-05 
i. 480-04 
2.152-03 
3.462-0 2 
5.820-01 
1.157+01 
2.981+02 
7.589+03 
1.389+05 
3.229+06 


EGAM 


1.545 

1.414 

1.299 

1.220 

1.186 

1.175 

1.172 

1.164 

1.145 


TCQOL 

(SEC) 

8.759-08 

1.001-07 

1.382-07 

2.697-07 

7.434-07 

2.793-06 

1.307-05 

5.675-05 

1.635-04 

5.246-04 


UNCLASSIFIED 


UNCLASSIFIED 


77 


DIAPHANOUS  DATA  CH2  GENERAL  ATONIC 


THETA (EV )  =  2.250 

SIGMA 

THETA  4  = 

2.630+13 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

E  ION 

ZBAR 

EGAN 

IGM/CM3)  IGM/CN3) 

(ERG/GM3) 

(ERGS/GM) 

( ERGS/GM J 

5.716-02  -2.862+00 

3.007+10 

1.186+12 

3.965+11 

0.139 

1.444 

5.746-03  -5.159+00 

3.389+09 

1.636+12 

7.510+11 

0.278 

1.361 

5.912-04  -7.433+00 

4.202+08 

2.511+12 

1.445+12 

0.539 

1.283 

7.764-05  -9.463+00 

6.530+07 

3.432+12 

2.170+12 

0.822 

1.245 

1.315-05  -1.124+01 

1.196+07 

3.978+12 

2.614+12 

0.970 

1.229 

2.406-06  -1.294+01 

2.289+06 

4.433+12 

3.006+12 

1.060 

1.215 

3.932-07  -1.475+01 

3.981+05 

5.207+12 

3.688+12 

1.193 

1.194 

6.148-08  -1.660+01 

6.520+04 

5.831+12 

4.240+12 

1.297 

1.182 

9.216-09  -1.850+01 

9.918+03 

6.052+12 

4.438+12 

1.331 

1.178 

1.171-09  -2.057+01 

1.277+03 

6.343+12 

4.707+12 

1.362 

1.172 

RHO 

LOG  RHO 

KROS 

KPLK 

LMUAKOS 

LMOAPLK 

TCOOL 

(  GM/CM3 } 

(GM/CM3) 

(CM2/GM) 

ICM2/GM) 

(  CM) 

(CM) 

(  SECS 

5.716-02 

-2.862+00 

9.405+04 

6.148+05 

1.860-04 

2.846-05 

1.834-08 

5.746-03 

-5.159+00 

6.778+04 

5.066+05 

2.568-03 

3.435-04 

3.070-08 

5.912-04 

-7.433+00 

4.062+04 

3.249+05 

4.164-02 

5.207-03 

7.348-08 

7.764-05 

-9.463+00 

1.489+04 

1.402+05 

8.647-01 

9.188-02 

2.328-07 

1.315-05 

-1.124+01 

4.089+03 

4.783+04 

1.860+01 

1.590+00 

7.908-07 

2.406-06 

-1.294+01 

1.111+03 

1.945+04 

3.742+02 

2.137+01 

2.167-06 

3.932-07 

-1.475+01 

2.571+02 

8.068+03 

9.893+03 

3.152+02 

6. 136-06 

6.148-08 

-1.660+01 

4.294+01 

2.462+03 

3.788+05 

6.608+03 

2.252-05 

9.216-09 

-1.850+01 

6.350+00 

9.219+02 

1.709+07 

1.177+05 

6.242-05 

1.171-09 

-2.057+01 

1.095+00 

6.079+02 

7.799+08 

1.405+06 

9.921-05 

UNCLASSIFIED 


UNCLASSIFIED 


78 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


THETA { EV )  =  3.400 

SIGMA 

THETA  4  = 

1.371+14 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

EION 

ZBAR 

EGAM 

{ GM/CM3 )  ( GM/CM3  > 

( ERG/CM3 ) 

(ERGS/GM) 

(ERGS/GM) 

3.421-02  -3.375+00 

3.420+10 

2.481+12 

9.809+11 

0.433 

1.403 

4.522-03  -5.399+00 

5.221+09 

3.537+12 

1.805+12 

0.655 

1.326 

6.660-04  -7.314+00 

8.780+08 

4.581+12 

2.603+12 

U.890 

1.288 

1.088-04  -9.126+00 

1.586+08 

5.536+12 

3.349+12 

1.089 

1.263 

1.910-05  -1.087+01 

2.986+07 

6.446+12 

4.102+12 

1.240 

1.242 

3.570-06  -1.254+01 

5.797+06 

7.068+12 

4.632+12 

1.328 

1.230 

6.174-07  -1.430+01 

1.039+06 

7.860+12 

5.337+12 

1.411 

1.214 

9.543-08  -1.616+01 

1.699+05 

9.302+12 

6.631+12 

1.552 

1.191 

1,380-08  -1.810+01 

2.553+04 

1.035+13 

7.574+12 

1.651 

1.179 

1.717-09  -2.018+01 

3.265+03 

1.130+13 

8.452+12 

1.725 

i  .168 

RHO 

LOG  RHO 

KROS 

KPLK 

LMOAROS 

LMDAPLK 

TCQQL 

( GM/CM3 ) 

(GM/CM3) 

(CM2/GM) 

{ CM2/GM) 

(CM) 

(CM) 

(SEC) 

3.421-02 

-3.375+00 

2.607+05 

5.231+05 

1.121-04 

5.588-05 

8.646-09 

4.522-03 

-5.399+00 

1.804+05 

3.062+05 

1.226-03 

7.222-04 

2.106-08 

6.660-04 

-7 -  314+00 

8.365+04 

1.505+05 

1.795-02 

9.978-03 

5.551-08 

1.088-04 

-9.126+00 

2.527+04 

7.340+04 

3.637-01 

1.252-01 

1.375-07 

1.910-05 

-1.087+01 

5.163+03 

3.584+04 

1.014+01 

1.460+00 

3.280-07 

3.570-06 

-1.254+01 

8.768+0 2 

2.034+04 

3.195+02 

1.377+01 

6.337-07 

6.174-07 

-i. 430+01 

1.830+02 

1.386+04 

8.849+03 

1.169+02 

1.034-06 

9.543-08 

-1.616+01 

3.448+01 

7.869+03 

3.039+05 

1.332+03 

2.155-06 

1.380-08 

-1 .810*01 

5.670+00 

4.224+03 

1.278+07 

1.715+04 

4.467-06 

i. 717-09 

-2. 018i-0i 

1. 078+00 

2.849+03 

5.402+08 

2.044+05 

7.234-06 

UNCLASSIFIED 


UNCLASSIFIED 


DIAPHANOUS  DATA 


79 

CH2  GENERAL  ATOMIC 

5.000  SIGMA  THETA  4  =  6.412*14 

LOG  RHO  PRESSURE  ENERGY  EION  Z8AR 

(GM/CM3J  (ERG/ CM3)  (ERGS/GH)  ( ERGS/GM ) 


THETA (EV) 
RHO 

( GM/CM3) 

3.566-02 
5. 6 80-03 
9.345-04 
1.600-04 
2.865-05 
5.238-06 
9.010-07 
1.412-07 
2.062-08 
2.647-09 


RHO 

IGM/CM3) 


-3.334+00 

-5.171+00 

-6.976+00 

-8.740+00 

-1.046+01 

-1.216+01 

-1.392+01 

-1.577+01 

-1.770+01 

-1.975+01 


LOG  RHO 
CGM/CM3) 


6.369+10 

1.125+10 

2.043+09 

3.810+08 

7.275+07 

1.405+07 

2.518+06 

4.159+05 

6.285+04 

8.135+03 


KROS 

CCM2/GM1 

1.855+05 
1.426+05 
6.948+04 
2.088+04 
4,331+03 
S. 000+02 
1.508+02 
2.803+01 
4.221+00 
7.221-01 


4.405+12 

5.786+12 

7.205+12 

8.467+12 

9.875+12 

1.134+13 

1.272+13 

1.478+13 

1.617+13 

1.654+13 


KPLK 
( CM2/GMI 

3.234+05 

1.819+05 

1.042+05 

7.259+04 

5.160+04 

3.437*04 

2.287+04 

1.056+04 

2.402+U3 

3.396+02 


1.726+12 

2.816*12 

3.926+12 

4.895+12 

6.U65+12 

7.314+12 

8.525+12 

1.033+13 

1.160+13 

1.193+13 


LMDARQS 
( CM ) 

1.511-04 

1.234-03 

1.540-02 

2.993-01 

8.060+00 

2.386+02 

7.357+03 

2.526+05 

1.149+07 

5.232*08 


0.741 

0.930 

1.131 

1.321 

1.475 

1.614 

1.724 

1.870 

1.970 

1.99b 


LMDAPLK 

(CM) 

8.670-05 
9.677-04 
1.027-02 
8.609-02 
6.765-01 
5.554+00 
4.852+Ui 
6.706+U2 
2.  19+04 
1,113+06 


EGAN 


1.405 
1.342 
1.303 
1.281 
1.257 
1.2  37 
1.220 
1.200 
1.188 
1.186 


TCOOL 
l  SEC) 

5.310-09 
1.240-08 
2.697-08 
4. 547-08 
7.461-08 
1.286-07 
2. 168-07 
5.459-07 
2.625-06 
1. 898-05 


3.566-02  -3.334+00 
5.680-03  -5.171+00 
9.345-04  -6.976+00 
1.600-04  -8.740+00 
2.865-05  -1.046+01 
5.238-06  -1.216+01 
9.010-07  -1.392+01 
1.412-07  -1.577+01 
2.062-08  -1.770+01 
2.647-09  -1.975*01 


UNCLASSIFIED 


UNCLASSIFIED 


80 


DIAPHANOUS  DATA 


CH2  GENERAL  ATOMIC 


IHLTA(bV)  =  7.000  SIGMA  THETA  4  =  2.463+15 


RHO 

LOG  RHO 

PRESSURE 

ENERGY 

EION 

ZBAR 

EGAM 

( GM/CM3 ) 

(GM/CM3) 

l ERG/CM3 I 

{ ERGS/GM1 

{ ERGS/GM ) 

4.490-02  -3.103+00 

1.273+11 

6.631+12 

2.376+12 

0.974 

1.428 

7.616-03  -4.878+00 

2.351+10 

8.547+12 

3.917+12 

1.149 

1 . 36 i. 

1.2bi-03  -6.660+UO 

4.354+09 

1.047+13 

5.370-12 

1.366 

1.325 

2.264-04  -8.393+00 

8.280+08 

1.244+13 

6.956+12 

1.546 

1.294 

4 .0  79-05  -1.011  +  01 

1.591+08 

1.441+13 

8.554+12 

1.716 

1.271 

7.490-06  -1.180+01 

3.087+07 

1.650+13 

1.032+13 

1.869 

1.25 

1.3IU-C6  -1.355+01 

5.578+06 

1.789+13 

1.150+13 

1.965 

1.238 

2.195-07  -1.533+01 

9.437+05 

1.832+13 

1.187+13 

1.993 

1.235 

3.367-08  -1. 721+01 

1.451+05 

1.842+13 

1.196+13 

1.999 

1.234 

4.3/6-09  -1.925+01 

1.885+U4 

1.844+13 

1.198+13 

2.000 

1.234 

RHO 

LOG  RHO 

KRUS 

KPLK 

LMDAROS 

LMDAPLK 

(COOL 

{ GM/CM3 ) 

(GM/CM3) 

(CM2/GM) 

(CM2/GM) 

CCMi 

(CM) 

(SEC) 

4.490-02 

-3.103+00 

9.396+04 

2.196+05 

2.370-04 

-.014-04 

3.064-09 

7.616-03 

-4.878+00 

8.351+04 

1.366+05 

1.572-03 

9.610-04 

6.348-09 

1.281-03 

-6.660+00 

4.855+04 

9.713+04 

1.608-02 

8.038-03 

1.094-08 

2.264-04 

-8.393+00 

1.438+04 

7.022+04 

3.072-01 

6.290-0 2 

798-08 

4.079-05 

-1.011+01 

3.139+03 

4.878+04 

7.810+00 

5. r-  25-01 

2.997-08 

7.490-06 

-1.180+01 

6.223+02 

2.342+04 

2.146+02 

5.702+UO 

7.152-08 

1.310-06 

-1.355+01 

8.469+01 

6.331+03 

9.015+03 

1.206+02 

2.868-07 

2.195-07 

-1.533+01 

1.058+01 

1.201+03 

4.306+05 

3.795+03 

1.549-06 

3.367-08 

-1.721+01 

1.614+00 

1.899+02 

1.839+07 

1.564+05 

9.845-06 

4.3/6-09 

-1. 925+01 

4.365-01 

2.487+01 

5.235+U8 

9.187+U6 

7.523-05 

UNCLASSIFIED 


UNCLASSIFIED 


81 

DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 

IHETAiEV)  =  10.000  SIGMA  THETA  4  =  1.026+16 


PRESSURE  ENERGY 
(ERG/CM3J  (ERGS/GM) 


RHO 

(GM/CM3) 


RHO 

t  GM/CM3 J 

6.208-02 

1.052-02 

1.884-03 

3.373-04 

6.234-05 

1.208-05 

2.2.2-06 

3.737-07 

5.747-08 

7.471-09 


LOG  RHO 
(GH/CM3 ) 


LOG  RHO 
(GM/CM3) 

-2.779+00 

-4.554+00 

-6.275+U0 

-7.994+00 

-9.683+00 

-1.132+01 

-1.303+01 

-1.480+01 

-1.667+01 

-1.871+01 


2.807+11 

5*225+10 

9.997+09 

1.919+09 

3.732+08 

7.384+07 

1.354+07 

2.300+06 

3.538+05 

4.591+04 


KROS 

{CM2/GMJ 

3.859+04 

3.494+04 

2.325+04 

8.266+03 

1.530+03 

2.109+02 

2.850+01 

4.366+00 

9.175-01 

3.629-01 


9.853+12 

1.261+13 

1.527+13 

1.790+13 

1.993+13 

2.081+13 

2.109+13 

2.118+13 

2.120+13 

2.121+13 


KPLK 

(CM2/GMJ 

1.488+05 

1.004+05 

6.613+0* 

3.721+04 

1.365+04 

3.523+03 

7.062+02 

1.237+02 

1.923+01 

2.511+00 


EION 

(ERGS/GM) 

3.070+12 

5.163+12 

7.304+12- 

9.371+12 

1.095+13 

1.163+13 

1.187+13 

1.194+13 

1.197+13 

1.198+13 


LMOAROS 

(CM) 

4.175-04 
2.720-03 
2.283-02 
3. 586-01 
1. 048+01 
3.925+02 
1.588+04 
6.129+05 
1.897+07 
3.689+08 


ZBAR 


1.230 

1.420 

1.587 

1.772 

1.917 

1.979 

1.996 

1.999 

2.000 

2.000 


LMDAPLK 

(CM) 

1.082-04 

9.467-04 

8.028-03 

7.967-02 

1.175+00 

2.350+Ul 

6.431+02 

2.164+^4 

9.046+05 

5.330+07 


EGAM 


1.459 

1.394 

1.348 

1.318 

1.300 

1.294 

1.291 

1.291 

1.290 

1.290 


TCOOL 

(SEC) 

1.613-09 

3.061-09 

5.625-09 

1.172- 08 
3.557-08 
1.439-07 
7.277-07 

4.172- 06 
2.686-05 
2.058-04 


6.208- 02  -2.779+00 
1.052-02  -4.554*00 
1.884-03  -6.275+00 
3.373-04  -7.994+00 
6.234-05  -9.683+00 

1.208- 05  -1.132+01 
2.202-06  -1.303+01 
3.737-07  -1.480+01 
5.747-08  -1.667+01 
7.471-09  -1.871+01 


UNCLASSIFIED 


UNCLASSIFIED 


02 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


THETA ( EV )  =  15.000 

SIGMA 

THETA  4  = 

5.194+16 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

E  ION 

ZBAR 

EGAM 

(GM/CM3)  (GM/CM3J 

( ERG/CM3 i 

(ERGS/GM) 

(ERGS/GM J 

9.543-02  -2.349+00 

7.161+11 

1.501+13 

3.751+12 

1.438 

1.500 

1.679-02  -4.087*00 

1.362+11 

1.883+13 

6.667+12 

1.635 

1.431 

3.053-03  -5.792+00 

2.629+10 

2.200+13 

9.079+12 

1.798 

1 . 392 

5.745-04  .462*00 

5.148+09 

2.413+13 

1.069+13 

1.911 

1.371 

1.117-04  -9.100*00 

1.020+09 

2.515+13 

1.145+13 

1.966 

1.363 

2.204-05  -1.072*01 

2.030+08 

2.558+13 

1.176+13 

1.992 

1.360 

4.040-06  -1.242+01 

3.728+07 

2.574+13 

1.189+13 

1.999 

1.359 

6.863-07  -1.419+01 

6.337+06 

2.580+13 

1.195+13 

2.000 

i  .358 

1.056-07  -1.606+01 

9.748+05 

2.582+13 

1.197+13 

2.000 

1.358 

1.372-08  -1.810+01 

1.267+05 

2.585+13 

1.200+13 

2.000 

1.357 

RHO 

LOG  RHO 

KRQS 

KPLK 

LMDAROS 

LMDAPLK 

TCOOL 

( GM/CM3 ) 

(GM/CM3J 

(CM2/GM) 

{ CM2/GM ) 

(CM) 

(CM) 

(SEC) 

9.543-02 

-2.349+00 

1.091+04 

8.170+04 

9.607-04 

1.283-04 

8.842-10 

1.679-02 

-4.087+00 

7.611+03 

4.632+04 

7.825-03 

1.286-03 

1.957-09 

3.053-03 

-5.792+00 

3.949+03 

2.260+04 

8*294-0 2 

1.449-02 

4.685-09 

5.745-04 

-7.462+00 

1.260+03 

7.548+03 

1.381+00 

2.306-01 

1.539-08 

1.117-04 

-9.100+00 

2.780+02 

1.906+03 

3.220+01 

4.697+00 

6.350-08 

2.204-05 

-1.072*01 

5.499+01 

4.331+02 

8.249+02 

1.047+02 

2.842-07 

4.040-06 

-1.242+01 

9.783+00 

8.388+01 

2.530*04 

2,951+03 

1.477-06 

6.863-07 

-1.419+01 

1.837+00 

1.492+01 

7.931+05 

9.762+04 

8.319-06 

1.056-07 

-1.606+01 

5.197-01 

2.721+00 

1.823+07 

3.481+06 

4.568-05 

1.372-08 

-1.810+01 

2.755-01 

7.739-01 

2.645+08 

9.416+07 

1.608-04 

UNCLASSIFIED 


UNCLASSIFIED 


83 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 

THETA (EV)  =  22.500  SIGMA  THETA  4  *  2.630+17 


RHO 

(GN/CM3) 

1.572-01 

2.833-02 

5.370-03 

1.039-03 

2.043-04 

4.043-05 

7.397-06 

1.240-06 

1.814-07 

2.205-08 


RHO 

( GM/CM3 ) 

1.572-01 

2.833-02 

5.370-03 

1.039-03 

2.043-04 

4.043-05 

7.397-06 

1.240-06 

1.814-07 

2.205-08 


LOG  RHO 
(GM/CM3) 

-1.850+00 

-3.564+00 

-5.227*00 

-6.870*00 

-8.496+00 

-1.012+01 

-1.181+01 

-1.360+01 

-1.552+01 

-1.763+01 


LOG  RHO 
(GM/CM3) 

-1.850*00 

-3.564+00 

-5.227*00 

-6.870+00 

-8.496*00 

-1.012+01 

-1.181+01 

-1.360+01 

-1.552+01 

-1.763+01 


PRESSURE 
( ERG/CM3 ) 

1.890+12 

3.636+11 

7.136+10 

1.411+10 

2.806+09 

5.592+08 

1.026+08 

1.737+07 

2.628+06 

3.346+05 


KROS 

(CM2/GM) 

3.246+03 
1.600+03 
5.328+0 2 
1.462+02 
3.357+01 
7.816+00 
1.974+00 
7.281-01 
3.792-01 
2.596-01 


ENERGY 

(ERGS/GM) 

2.174+13 

2.644+13 

2.954+13 

3.126+13 

3.210+13 

3.260+13 

3.320+13 

3.560+13 

4.487+13 

5.794+13 


KPLK 

(CN2/GH) 

3.662+04 

1.626+04 

5.382+03 

1.527+03 

4.199+02 

1.632+02 

1.036+02 

8.544+01 

5.661+01 

1.963+01 


EION 

(ERGS/GM) 

3.709+12 

7.188+12 

9.612+12 

1.088+13 

1.150+13 

1.185+13 

1.238+13 

1.460+13 

2,314+13 

3.517+13 


LMDAROS 

(CM) 

1.959-03 

2.206-02 

3.495-01 

6.586+00 

1.458+02 

3.164*03 

6.847+04 

1.107+06 

1.454*07 

1.747+08 


ZBAR 


1.604 

1.780 

1.878 

1.942 

1.975 

1.996 

2.005 

2.033 

2.139 

2.287 


LMDAPLK 

(CM) 

1.737-04 

2.171-03 

3.460-02 

6.305-01 

1.166+01 

1.515+02 

1.305+03 

9.436+03 

9.739+04 

2.311+06 


EGAM 


1.553 

1.485 

1.450 

1.435 

1.428 

1.424 

1.418 

1.393 

1.323 

1.262 


TCOOL 
( SEC  ) 

5.645-10 
1.546-09 
5.219-09 
1.946-08 
7.269-08 
1.899-07 
3.046-07 
3. 962-07 
7.536-07 
2. 806-06 


UNCLASSIFIED 


UNCLASSIFIED 


Qh 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


THETA (  EV I  =  34.000 

SIGMA 

THETA  4  « 

1.371+18 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

EION 

ZBAR 

EGAM 

( GM/CM3 )  CGM/CM3 J 

{ ERG/CM3) 

(ERGS/GM) 

(ERGS/GM) 

- 

2.731-01  -1.298+00 

5.173+12 

3.164+13 

3.230+12 

1.715 

1.599 

5.035-02  -2.989+00 

1.005+12 

3.723+13 

7.291+12 

1.860 

1.536 

9.700-03  -4.536+00 

1.984+11 

4.063+13 

9.948+12 

1.931 

1.503 

1.892-03  -6.270*00 

3.934+10 

4.356+13 

1.236+13 

1.981 

1.477 

3.651-04  -7.915+00 

7.775+09 

4.916+13 

1.721+13 

2.053 

1.433 

6.857-05  -9.588+00 

1.524+09 

6.Q65+13 

2.731+13 

2.186 

1.366 

RHO 

LOG  RHO 

KRQS 

KPLK 

LMDAROS 

LMDAPLK 

TCOOL 

(GM/CM3) 

(GM/CM3) 

(CM2/GM) 

(CM2/GM) 

(CM) 

(CM) 

(SEC) 

2.731-01 

-1.298+00 

1.587+03 

1 . 523+<.4 

2.308-03 

2.404-04 

3.788-10 

5.035-02 

-2. 989+00 

5.617+02 

6.282+03 

3.535-02 

3.161-03 

1.080-09 

9.7C0-03 

-4.636+00 

1.722+0 2 

3.116+03 

5.986-01 

3.308-02 

2.377-L9 

1.892-03 

-6.270+00 

4.596+01 

2.079+03 

1.150*01 

2.543-01 

3.820-09 

3.651-04 

-7.915+00 

1.441+01 

1.617+03 

1.901+02 

1.694+00 

5.543-09 

6.857-05 

-9.588+00 

5.226+00 

1.029+03 

2.790+03 

1.417+01 

1.075-08 

UNCLASSIFIED 


UNCLASSIFIED 


05 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


\ 

THETA! EV)  =  50.000 

SIGMA 

THETA  4  * 

6.412+18 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

E!ON 

ZBAR 

EGAM 

IGM/CM3 )  (GM/CM3) 

(ERG/ CM3 I 

(ERGS/GM) 

(ERGS/GM) 

4.678-01  -7.597-01 

1.337*13 

4.608+13 

3.215+12 

1.785 

1.620 

8.651-02  -2.447*00 

2.601+12 

5.495+13 

9.844+12 

1.931 

1.547 

1.616-02  -4.125*00 

5.086+11 

6.582+13 

1.861+13 

2.068 

1.478 

3.004-03  -5.808*00 

9.937+10 

8.089+13 

3.126+13 

2.224 

1.409 

) 

5.616-04  -7.485*00 

1.947*10 

9.689+13 

4.488+13 

2.380 

1.358 

1.055-04  -9.157*00 

3.825*09 

1.140+14 

5.964+13 

2.533 

1.318 

*■ 

RHO  LOG  RHO 

KROS 

KPLK 

LMDAROS 

LMDAPLK 

TCOOL 

! GM/CM3 )  (GM/CM3) 

CCM2/GMJ 

(CN2/GM) 

(CM) 

(CM) 

(SEC) 

4.678-01  -7.597-01 

1.694*03 

1.224+04 

1.262-03 

1.747-04 

1.468-10 

y 

8.651-02  -2.447*00 

5.590+02 

8.448+03 

2.068-02 

1.368-03 

2.536-10 

1.616-02  -4,125*00 

1.932+02 

6.414*03 

3.203-01 

9.649-03 

4.001-10 

3.004-03  -5.808*00 

6.124+01 

4.031+03 

5.436+00 

8.258-02 

7.823-10 

5.616-04  -7.485+00 

1.751+01 

2.201+03 

1.017+02 

8.089-01 

1.716-09 

1.055-04  -9.157+00 

4.801+00 

9.266+02 

1.974+03 

1.023+01 

4.797-09 

UNCLASSIFIED 


UNCLASSIFIED 


86 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


THETA (EV)  =  70.000 

SIGMA 

THETA  4  = 

2.463+19 

-  -  -  -  - 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

E  ION 

ZBAR 

EGAM 

( GM/CM3 )  (GM/CM3) 

(ERG/ CM3) 

(ERGS/GM) 

C ERGS/GM) 

7.098-01  -3.427-01 

3.007+13 

7.054+13 

6.991+12 

1 «  *'  49 

1.601 

1.303-01  -2. 038+00 

5.846+12 

8.921+13 

2.191+13 

2.124 

1.503 

2.402-02  -3.729+00 

1.140+12 

1.091+14 

3.789+13 

2.304 

1.435 

4.465-03  -5.411+00 

2.231+11 

1.297+14 

5.472+13 

2.479 

1.385 

8.508-04  -7.069+00 

4.402+10 

1.438+14 

6.622+13 

2.602 

1.360 

1.674-04  -8.695+00 

8.764+09 

1.491+14 

7.055+13 

2.644 

1.351 

RHO  LOG  RHO 

KRQS 

KPLK 

LMOAROS 

LMDAPLK 

TCOOL 

( GM/CH3 )  (GM/CM3) 

iCH2/GH) 

ICM2/GM) 

(CM ) 

(CM) 

(SEC) 

7.098-01  -3.427-01 

2.798+03 

1.354+04 

5.036-04 

1.040-04 

5.286-11 

1.303-01  -2.038+00 

1.026+03 

9.541+03 

7.483-03 

8.043-04 

9.489-11 

2.402-02  -3.729+00 

3.027+02 

5.824+C3 

1.375-01 

7.148-03 

1.901-10 

4.465-03  -5.411+00 

7.233+01 

2.723+03 

3.097+00 

8.224-02 

4.833-10 

8.508-04  -7.069+00 

1.456+01 

8.489+02 

8.072+01 

1.384+00 

1.719-09 

1.674-04  -8.695+00 

2.771+00 

1.959+02 

2.156+03 

3.048+01 

7.721-09 

UNCLASSIFIED 


UNCLASSIFIED 


8? 


DIAPHANOUS  DATA 


CK2  GENERAL  ATOMIC 


THETA (EV) 
RHO 

(GM/CM3  ) 

1.105+00 

2.035-01 

3.760-02 

7.243-03 

1.430-03 

2.846-04 


RHO 

(GM/CM3) 

1. 105+00 
2.035-01 
3.760-02 
7.243-03 
1.430-03 
2.846-04 


=  103.000 

LOG  RHO 
(GM/CM3) 

9.978-02 
-1.592+00 
-3.281+00 
-4. 928+00 
-6.550+00 
-8.164+00 


LOG  RHO 
(GM/CM3) 

9.978-02 
-1.592+00 
-3.281+00 
-4.928+00 
-6.550+00 
—  8  « 164+00 


SIGMA 

PRESSURE 
t ERG/ CM3 I 

7.115+13 

1.387+13 

2.711+12 

5.364+11 

1.069+11 

2.135+10 


KROS 

(CM2/GH) 

3.347+03 

1.340+03 

3.361+02 

5.853+01 

9.215+00 

1.848+00 


THETA  4  = 

ENERGY 

(ERGS/GM) 

1.129+14 

1.397+14 

1.638+14 

1.773+14 

1.822+14 

1.840+14 


KPLK 

(CM2/GMI 

1.071+04 

6.121+03 

2.602+03 

7.231+02 

1.655+02 

3.500+01 


1.026+20 
El  ON 

(ERGS/GM) 

1.631+13 

3.742+13 

5.563+13 

6.618+13 

7.007413 

7.143+13 


LHDAROS 

(CM) 

2.704-04 

3.668-03 

7.913-02 

2.359+00 

7.587+01 

1.901+03 


ZBAR 


2.138 
2.322 
2.513 
2.609 
2 . 643 
2.656 


LMDAPLK 

(C.M) 

8.449t05 

8.030-04 

1.022-02 

1.909-01 

4.224+00 

1.004+02 


EGAM 


1.570 

1.488 

1.440 

1.418 

1.410 

1.408 


TCOOL 

(SEC) 

2.568-11 

5.561-li 

1.534-10 

5.974-10 

2.682-09 

1.281-08 


UNCLASSIFIED 


r-4  nO  -4  CsJ  LTV 


UNCLASSIFIED 

88 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


THETA (EV)  =  150.000 

SIGMA 

THETA  4  = 

5.194+20 

RHO  LOG  RHO 

PRESSURE 

ENERGY 

EION 

2BAR 

EGAM 

(GM/CM3)  IGM/CM3) 

( ERG/CM3 ) 

i ERGS/GM) 

I ERGS/GM J 

1.855+00  6.179-01 

1.907+14 

1.812+14 

2.705+13 

2.340 

1.5&7 

3.476-01  -1.057+00 

3.742+13 

2.112+14 

4.975+13 

2.497 

1.510 

6.681-02  -2.706+00 

7.400+12 

2.284+14 

6.223+13 

2  a  599 

1*485 

1.319-02  -4.329+00 

1.475+12 

2.357+14 

6.796+13 

2.633 

1.474 

2.618-03  -5.945+00 

2.943+11 

2.390+14 

7.039+13 

2.653 

1.470 

5.223-04  -7.557+00 

5.882+10 

2.404+14 

7.145+13 

2.659 

1.468 

RHQ 

LOG  RHO 

KROS 

KPLK 

LHOAROS 

LMDAPLK 

TCQOL 

(GM/CM3) 

( GM/CM3 I 

(CM2/GM) 

( CM2/GM) 

(CM) 

(CM) 

(SEC) 

1.855+00 

6.179-01 

1 -Q2 1*03 

5.042+03 

1*069^04 

i 1 i 

3. 476-01 

-1.057+00 

3.924+02 

1.960+03 

7.332-03 

1.468-03 

5.187-11 

6.681-02 

-2.706+00 

1.078+02 

5.811+02 

1. 368-01 

2.576-02 

1.892-10 

1.319-02 

-4.329+00 

2.308+01 

1.375+U2 

3.286+00 

5.517-01 

8. 253-10 

2.61 8—03 

rr  rw  •  i~  .  ^ 

T*  f'fT'fOO 

2.982+01 

8.051+01 

1.281+01 

3.859-09 

5.223-04 

-7.557+00 

1.208+00 

6.165+00 

1.584+03 

3.105+02 

1.877-08 

UNCLASSIFIED 


09 


diaphanous  data  cm2  general  atumc 


IHETAIEV) 

-  225.000 

SIGMA 

THETA  4  = 

2.630+21 

RHO 

LOG  RHO 

PRESSURE 

ENERGY 

E  ION 

ZBAR 

EGAH 

(GM/CM3) 

( GM/CM3 ) 

(ERG/CM3) 

(ERGS/GM) 

l ERGS/GM) 

3.220+00 

1.170+00 

5.168+14 

2.725+14 

3.180+13 

2.476 

1.589 

6.180-01 

-4. 813-01 

1.022+14 

2.995+14 

5.154+13 

2.580 

1.552 

1.219-01 

-2.105+00 

2.035+13 

3.132+14 

6.270+13 

2.617 

1.533 

2.410-02 

-3.725+00 

4.058+12 

3.204+14 

6.784+13 

2.647 

1.525 

4.807-03 

-5.336+00 

8.109+11 

3.233+14 

7.027+13 

2.654 

1.522 

9.589-04 

-6.950+00 

1.621+11 

3.249+14 

7.138+13 

2.661 

1.520 

RHO 

LOG  RHO 

KROS 

KPLK 

LMDAROS 

LMDAPLK 

TCOOL 

( GM/CM3 ) 

( GM/CM3 ) 

ICM2/GM) 

l CM2/GM) 

(CM) 

(CM) 

(  SEC ) 

3.220+00 

1.170+00 

1.738+02 

1.657+03 

1.786-03 

1.874-04 

1.563-11 

6.180-01 

-4.813-01 

6.450+01 

5.769+02 

2.509-02 

2.805-03 

4.936-11 

1.219-01 

-2.105+00 

1. 553+01 

1.350+02 

5.283-01 

6.0/6-02 

2.205-10 

2.410-02 

-3.725+00 

4.029+00 

3.102+01 

1.030+01 

1 .337+00 

9. 819-10 

4.807-03 

-5. 338+00 

1.182+00 

6.439+00 

1.759+02 

3.230+01 

4.774-09 

9.589-04 

-6.950+00 

4.939-01 

1.323+00 

2.111+03 

7.883+02 

2.335-08 

UNCLASSIFIED 


UNCLASSIFIED 


90 


DIAPHANOUS  data  CH2  GENERAL  ATOMIC 


THETAtEW) 

=  340.000 

SIGMA 

THETA  4  = 

1.371+2 2 

RHO 

LOG  RHO 

PRESSURE 

EMERGV 

eion 

?8AR 

EGAM 

IGM/CM3J 

( GM/CM3 ! 

IERG/CM3) 

! ERGS/GM) 

{ ERGS/GM J 

5.731+00 

1.746+00 

1.433+15 

4.038*14 

2.872*13 

2.584 

1.619 

1.132+00 

1.242-01 

2.857+14 

4.290*14 

5.053+13 

2.616 

1.588 

2.244-01 

-1.494+00 

5.699+13 

4.428+14 

6.185+13 

2.640 

1.574 

4.474-02 

-3.107+00 

1.139+13 

4.492+14 

6.740+13 

2.648 

1.567 

8.918-03 

-4.720+00 

2.275+12 

4.528*14 

7.005+13 

2.657 

1.564 

1.781-03 

-6.330+00 

4.549+11 

4.544+14 

7.128+13 

2.661 

1.562 

RHO 

LOG  RHO 

KROS 

KPLK 

LMOAROS 

LMDAPLK 

TCOOL 

(GM/CM3* 

IGM/CM3) 

(CM2/GM) 

ICM2/GM1 

ICM1 

(CM) 

(SECI 

5.731+00 

1.746+00 

4.031+01 

3.803+02 

4.328-03 

4.587-04 

1.936-11 

1.132+00 

1.242-01 

1.149+01 

1,561+02 

7.688-02 

5.656-03 

5.009-11 

2.244-01 

-1.494*00 

3.116+00 

3.593+01 

1.430+00 

1.240-C1 

2.247-10 

4.474-02 

-3.107+00 

1.010+00 

7,587+00 

2. 213+01 

2.946+00 

1.080-09 

8.918-03 

-4.720+00 

4.688-01 

1.581+00 

2.392+02 

7.094+Ul 

5o 223-09 

1.781-03 

-6.330+00 

3.033-01 

3.218-01 

1.851+03 

1.744*03 

2.575-08 

UNCLASSIFIED 


Ui 


UNCLASSIFIED 


91 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 

THETA <EV)  =  500.000  SIGKA  THETA  4  =  6.412+22 

RHO  LOG  RHO  PRESSURE  ENERGY  E ION  ZBAR  EGAM 

CGM/CM3)  CGM/CM3)  ( ERG/CM3  )  I ERGS/GM)  tCRGS/GM J 


1.013+01  2.315+00 

2.012+00  6.992-01 

3.997-01  -9.169-01 
7.952-0 2  -2.532+00 
.590-02  -4.142+00 
.175-03  -5.752+00 


3.749+15 

7.484+14 

1.494+14 

2.984+13 

5.967*12 

1.193+12 


5.796+14 
6.061+14 
6.214+14 
6.298+14 
6.329+14 
6 • 348+14 


2.439+13 
4.820+13 
6. 073+13 
6.688+13 
6.981+13 
7.117*13 


2.608 

2.625 

2.643 

2.657 

2.659 

2.662 


1.639 

1.614 

1.601 

1.596 

1.593 

1.592 


RHO 

( GM/CM3 ) 

1.013+01 

2.012+00 

3.997-01 

7.952-02 

1.590-02 

3.175-03 


LOG  RHO 
(GM/CM3) 

2.315*00 

6.992-01 

-9.169-01 

-2.532+00 

-4.142+00 

-5.752+00 


KRQS 
(CM2/GM } 

1.092+01 
3. 196+00 
1.068+00 
4.887-01 
3.117-01 
2.532-01 


KPLK 

(CM2/GH) 

i. 996+02 
4.964+Ul 
1.103+01 
2.320+00 
4.735-U1 
9.614-U2 


LMOAROS 

(CM) 

9.039-03 
i. 555-01 
2.341+00 
2.573+01 
2.018+02 
1.244+03 


LMDAPLK 

(CM) 

4.946-04 

1.001-02 

2.269-01 

5.420+00 

1.328+02 

3.276+03 


TCOOL 

(SEC) 

1.132-11 

4.761-11 

2.197-10 

1.058-09 

5.210-09 

2.574-08 


UNCLASSIFIED 


UNCLASSIFIED 


92 


DIAPHANOUS  DATA  CH2  GENERAL  ATOMIC 


IHETA(EV) 

=  700.000 

SIGMA 

THETA  4  == 

2.463+23 

RHO 

LOG  RHO 

PRESSURE 

ENERGY 

El  ON 

ZBAR 

EGAM 

1 GM/CM3 ) 

(GM/CM3) 

{  ERG/CM3  3 

(ERGS/GM} 

(ERGS/GM) 

I .67 1  +  01 

2.816+00 

8.685+15 

7.985+14 

1.888+13 

2.618 

1.651 

3.326+00 

1  ,  202+uO 

1.735+15 

8*277+14 

4.570+13 

2.629 

1.630 

6.613-01 

-4. 128-01 

3,464+14 

8.448+14 

5.957+13 

2.644 

1.620 

1.317-01 

-2.D27+U0 

6.919+13 

8.545+14 

6.635+13 

2.658 

1.615 

RHO 

LOG  RHO 

KROS 

KPLK 

LMDAROS 

LMOAPLK 

TCOOL 

( GM/CM3 ) 

(GM/CM3) 

(CM2/GM! 

(CM2/GM) 

(CM) 

(CM) 

(  SEC) 

1.6/1+01 

2.816+00 

4.208+00 

8.175+01 

1.422-02 

7.321-04 

9.913-12 

3.323+00 

1.202+00 

1.383+00 

1.937+01 

2.172-01 

1.551-02 

4.337-11 

6,618-01 

-4.128-01 

5.795-01 

4.139+00 

2,607+00 

3.651-01 

2.071-10 

3 .317-01 

-2.027+00 

3.439-U1 

n  Ten _ Al 

O ■ I>U”Vi 

O  o.i  Q  -L  A  T 

8.678+00 

9,910-1 

UNCLASSIFIED 


UNCLASSIFIED 


DIAPHANOUS  DA f A 
IriLIA(EV)  =  lUCO.dOO 


93 


CH2  GENERAL  ATOMIC 

SIGMA  THETA  4  =  I. 026+24 

LOG  RHO  PRESSURE  ENERGY  El  ON  ZBAR 

CGM/CM3)  l ERG/GM3J  1 6RGS/GM )  (ERGS/GM) 


RHO 

IGH/CM3 ) 

2.846+01 
5 .676+00 
1.130+00 
2./ 49-01 


RHO 

(GM/CM3  ) 

2.846+01 

6.676+00 

1.13C+00 

2.249-01 


3.348+00 

1.736+00 

1.220-01 

-1.492+00 


LOG  RIG 
(o M/CM  3  ) 

3.348+00 

1.736*00 

1.220-Ul 

-1.492+00 


2.117+16 
4.230+15 
2.450+14 
1 • 688+ 14 


KROS 
t  C  M2/GM ) 

1.862+00 

7.231-01 

3.904-01 

2.812-01 


1.128+15 

1.161+15 

1.180+15 

1.192*15 


KPLK 
l CM2/GM) 

3.309+01 

7.361+00 

1.579+00 

3.307-Ul 


1.226+13 

4.268+13 

5.817+13 

6.574+13 


LMDAROS 
(CM  J 

1.887-02 

2.437-01 

2.268+00 

1.582+01 


2.625 

2.633 

2.645 

2.658 


lmdaplk 

t  CM ) 

1.062-03 

2.394-02 

5.606-01 

1.345+01 


EGAM 


1.659 

1.642 

1.634 

1.630 


TCQOL 
(  SEC) 

8. 307-12 
3.843-11 
1 • 821-1 Q 
8.781-10 


UNCLASSIFIED 


UNCLASSIFIED 


9h 


DIAPHANOUS  DATA 


CH2  GENERAL  ATOMIC 

=  2250* 000  SIGMA  THETA  A  *  2.630*25 


ZBAft  EGAN 


THETA  t  tV ) 

RHO 

(GM/CM3) 

9.580*01 

1.897*01 

3.796*00 

7.588-01 


RHO 

( GM/CM3 ) 

9.580+01 

1.897*01 

j.?96»0:. 

7.588-01 


LOG  RHO 
IGM/CM3) 

A. 562+00 
2.943+00 
1.334*00 
2.760-01 


LOG  RHO 
(GM/CM3J 

4.562+00 
2.943+00 
1.334  +  00 
2.760-01 


PRESSURE 
l  ERG/ CM3 1 

1,606+17 

3.204+16 

6.409+15 

1,282+15 


KROS 

CCH2/GM) 

5.619-01 

3.378-Ox 

2.633-01 

2.300-01 


ENERGY 
{ tRGS/GM) 

2.509+15 

2.568+15 

2.587+15 

2.598+15 


KPLK 
I CH2/GM1 

4.843+00 

1.019+00 

2.139-01 

4.347-02 


E  ION 

(ERGS/GM) 

6.202+12 

3.441*13 

5.452*13 

6.416+13 


lhdaros 

(CM) 

1.058-02 

1,561-01 

L.GCO+OO 

5.537+011 


2.632 

2.659 

2.657 

2.658 


LMDAPLK 

(CM) 

2.156-03 
5.175-02 
1. 232+OC 
3. 031+u 1 


1.660 

1.658 

1.653 

1.650 


TCOOL 

(SEC) 

4.926-12 
2. 397-11 
1.150-10 
b. 681-10 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


